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1 . 0 SUMMARY 


REPKODiJCin:: :: 

ORIOINA- 


• the 
is i?OOR 


The objectives of the Blade Row Dynamic Digital Compressor Program were 
twofold. Firstly, the General Electric developed pitch-line, blade row, time 
marching, digital compression component stability model was adapted to the 
J85-13 engine configuration. In particular, the J85-13 compressor, including 
the combustor volume to the turbine diaphragm, was modeled taking into account 
the variable IGV geometry and the third, fourth, and fifth stage bleeds. The 
clean inlet performance of the compressor component was reproduced for two 
engines including dynamic indication of the surge line. This prediction for 
surge was accomplished by developing a stability criterion based upon the 
derivative of flow rate within the blade rows as compared to the derivative 
of the flow imposed by the throttling process boundary condition. Secondly, 
the clean inlet flow compressor model was modified to a parallel compressor 
configuration to permit imposing total-pressure, total-temperature, and com- 
bined total-pressure and total-temperature distorted upstream boundary con- 
ditions. The flow split between the sectors was determined by a method which 
simultaneously satisfied the inlet boundary conditions and the parallel- 
compressor uniform-static-pressure assumption imposed at the entrance to the 
combustor volume. It was anticipated that static-pressure gradients should 
be minimal at this location due to the low Mach number of the flow. Since 
the input to this model is the clean-inlet-flow blade row characteristics 
(relative total-pressure loss coefficients and deviation angles) as functions 
of incidence angle and corrected speed, it was necessary to develop a pro- 
cedure for defining blade row characteristics at corrected speeds other than 
those for which clean- inlet-flow data existed. In this manner good simu- 
lations with total- temperature distorted inlet flow boundary conditions could 
be obtained. This procedure makes use of spline curve fits of the blade char- 
acteristic data as a function of corrected speed within the data range. This 
permits interpolations to be carried out inside the range of data, while out- 
side the data range, the blade characteristics are obtained from linear 
extrapolations of the spline at its end points. 

Simulations of 25 180-degree 1/rev circumferential distortions were 
carried out for the following types of distortion patterns and the indicated 
speed range: 

A. Pure Total Pressure (80%-100% N//^) 

B. Pure Total Temperature (87%-100% N//¥) 

C. Combined Total Pressure and Total Temperature 

180° Opposed (87%-100% N//e) 

180° Coincident (87%-100% N//0) 

90° Overlapped (87%-100% N//0^) 

The loss in surge pressure ratio trends were correctly predicted at all 
speeds. However, the model exhibits characteristics similar to other parallel 
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compressor models, that is, it predicts the loss in surge pressure ratio 
accurately when the compressor speed line is near vertical and over predicts 
the loss in surge pressure ratio when the compressor speed line has a low 
slope. This deficiency is being studied and atteaqpts to rectify It, while 
still maintaining the major concepts of the classical parallel-compressor 
model, are being carried out in a continuing study. 
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2.0 INTRODUCTION 


Engine inlet total-pressure and total- temperature distortions lead to 
well-documented losses of surge pressure ratio in the compression systems of 
turbofan and turbojet engines. One approach to obtaining a more detailed 
understanding of the internal flow mechanisms which cause this loss in surge 
pressure ratio is to develop computer simulations of the compression systems 
that will not only permit studying the effects of steady-state spatial dis- 
tortions, but provide the means for rapidly and efficiently screening the 
effects of many types of distortion. Such simulations are also useful for 
determining the potential effects of design modifications on compression 
system stability. 

However, if such models are to achieve their ultimate capability as a 
design and evaluation tool, it is necessary to develop and validate the 
capability of the models against existing bodies of test data. This step 
establishes confidence in the model, helps to define its range of validity, 
and insures that it is understood. Toward this goal, this report presents 
the results of an analytical investigation in which an existing generalized 
compression component computer model was modified to simulate the performance 
characteristics of the J85-13 compressor system operating without inlet dis- 
tortion. The model was then modified to predict the effects of inlet dis- 
tortion on the stability and internal flow characteristics of the compressor 
system. The resulting predictions were compared to the measured effects. 

An existing one-dimensional pitch-line, dynamic digital model for com- 
pression components, developed by the Aircraft Engine Group of the General 
Electric Company, was used as the basis for the investigation. Because the 
breakdown of flow in a compression system is an inherently unsteady aero- 
dynamic phenomena which typically manifests itself as a rotating stall or 
surge, it was postulated that a time- dependent model would offer a unique 
approach for studying the factors affecting compression system stability. 
Certainly, one could question the value of a one-dimensional pitch-line 
model for investigating stability problems since it is known that rotating 
stall and surge are multi- dimensional events which generally initiate and 
propagate in the hub or tip regions of the blading. However, for aerodynamic 
stability studies it is not necessary to be able to detail the propagation 
velocity, number of cells, size, etc., of a rotating stall or the spatial 
distribution of the wave front of the surge pulse as long as the conditions 
under which the aerodynamic instability would occur can be determined. A 
properly constructed one-dimensional model should have the ability to pre- 
dict the circumstances under which a disturbance will change from being 
attenuated to being propagated. Examination of a stability criterion derived 
by Jansen (Reference 1) supports this contention. 

Another way to convey the same information is to examine the eigenvalues 
of the characteristic equation derived from the Jacobian matrix of equations 
describing the aerodynamic performance of the compression system. This 
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approach is currently under Investigation in a parallel study (Contract No. 
NAS3-19854). In it, a change from negative to positive real eigenvalues 
will lead to the propagation of disturbances - a problem akin to the problem 
of acoustic propagation in wave guides. The use of this method has been 
discussed by Danlele, Blaha, and Seldner (Reference 2). 

Confidence in the basic approach used in the present study had been 
established in a study documented by Ruegg (Reference 3). In it, the propa- 
gation velocity of waves in a duct were studied using the basic equations of 
the model and the calculatlonal technique. The results were compared with 
method of characteristics results. Analysis showed that the model produced 
accurate predictions of wave speed, thermodynamic properties, and flow 
parameters. Resonances in ducts were also studied and the resonant frequen- 
cies compared well with those predicted by acoustic theory. Proper location 
of nodal points with frequency was noticed as well as proper qualitative 
changes in wave amplitude with frequency. 

In another study, both a two— stage fan and a nine— stage compressor were 
simulated with clean inlet conditions. The clean inlet maps were accurately 
reproduced and the surge lines predicted. Time dependent inlet and exit 
boundary conditions were Imposed without creating numerical instabilities, 
and the results Indicated proper qualitative response. 

With this background, it was felt that the generalized model which had 
demonstrated the ability to properly calculate the state of the fluid in 
ducts and blade rows, was in a sufficient state of development to adapt it 
to the "real world" stability problems encountered in compression systems. 

The compressor model was divided into volumes, one blade row per volume, 
except for free volumes whose lengths were chosen to be commensurate with 
the longest axial length blade row. The equations of change (conservation 
of mass, momentum, and energy) integrated once over the volume to give the 
macrobalances were written in a form that permitted the determination of 
time derivatives of density, physical flow, and the product of density times 
entropy as fimctions of space variables. Time-dependent solution was accom- 
plished by substituting the time derivatives in a Taylor series to give an 
estimate of the three above mentioned state variables at the next increment 
in time. This was continued until the solution had settled out if the bound- 
ary conditions were not time dependent, or until stationary behavior was 
reached if periodic boundary conditions were Imposed, or for any portion of 
a transient for boundary conditions not previously stated. 

It is important to note that the momentum and energy macrobalance equa- 
tions from which the time derivatives of flow and density-entropy product are 
determined contain a blade force term and an entropy production term, respec- 
tively. The blade force was obtained from resolving the forces that act on a 
blade including the tangential force obtained from the Euler turbine equation, 
and the blade drag force which was related to the losses within a blade row. 
The entropy production term was also related to the losses that develop with- 
in a blade row. In this generalized model, it is the presence of these terms 
which determine that a volume is treated as a blade row volume. The absence 
of the blade force and entropy production terms Indicates a free volume 
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indicative of lossless duct flow. Hence, it was through the blade force and 
entropy production terms that the clean inlet performance of the compressor 
was input. The input of the performance was accomplished through the rela- 
tive total-pressure loss coefficient and the deviation angle as a function of 
incidence angle at each corrected speed for each rotor. The stators were as- 
sumed to be lossless although a constant deviation angle was assigned to them 
The relative total-pressure loss coefficient and the rotor deviation angle 
were determined from the steady-state stage characteristics. Because the 
method completely specifies the flow condition, program output Includes the 
total temperatures and total pressures of individual stages, stage character- 
istics, velocity diagram information at both the inlet and exit of each blade 
row, and diffusion factors. 

The present investigation was divided into two parts - (1) Clean Inlet 
Model and (2) Distorted Inlet Model. For the clean inlet model, the existing 
General Electric Dynamic Digital Blade Row Compression Component Stability 
Model was modified to represent the NASA-Lewis Research Center J85-13 com- 
pressor configurations and stage performance characteristics for undistorted 
inlet flow conditions. The resulting model accounted for the scheduled 
changes of variable IGV angle and variable third, fourth, and fifth stage 
bleed flow with changes in corrected speed. The performance characteristics 
of each blade row were determined from stage stacking procedures and expressed 
in terms of Incidence angle, deviation angles, and loss coefficient variations 
or each blade row. These values were determined from conq>re8sor Interstage 
data furnished by the Lewis Research Center which consisted of a hub, mean 
and tip radius total-temperature, total-pressure and tip wall-static-pressure 
measurements for the pressure ratio and flow ranges of each speed line. The 
resulting model also incorporated an Improved stability criterion, based on 
the self developing unsteady internal flows generated near surge. 

Two clean inlet compressor maps were generated using the computer model 
The maps consisted of four speed lines (80, 87, 94, and 100 percent corrected 
speed) for the first, or "Moss" engine (Reference 4), three speed lines (87 
94 and 100 percent corrected speed) for the second, or "Mehalic" engine 
(Reference 5), and the predicted surge points for each speed line. Verifies- 
tlon of the computer model was made based on the comparisons achieved between 
the predicted maps and the corresponding experimental maps presented in the 
references. These comparisons were based on the accuracy of predicting the 
pressure ratio flow characteristics, and the surge point for each speed 
investigated. Following this verification, the predicted stage velocity 
diagrams, blade and stator incidence angles, diffusion factors, and loss 
coefficients at surge and at two points below surge for each corrected speed 
line were analytically determined. 

The computer model was then modified to a multi-sector parallel compressor 
configuration to accept circumferentially distorted inlet flow conditions. 

This modification was accomplished by dividing the inlet annulus into sectors 
the number of which depended on the manner in which the state properties varied 
circumferentially at the compressor inlet. The maximum number of sectors was 
limited to 12, and the minimum number was limited to two. The variable IGV 
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and bleed schedules, and the stage performance characteristics used for each 
distorted sector were the same as those used for the corresponding compressor 
with undistorted inlet flow. As the rotor blades moved from one sector to 
the next it was assumed that both the rotor inlet and exit velocity diagrams 
shifted instantaneously to a new equilibrium condition representative of the 
next sector. No account was made of the possible effects of the rotor blade 
unsteady response characteristics. This approach was deemed appropriate for 
the J85 compressor operating with 180 degree extent distortion, where the 
broad extent of the distorted sector coupled with the short chordedrot^ 
blades resulted in low values of rotor reduced frequency (w < 0.17). The 
analysis of Schorr and Heddy (Reference 6) indicates the essentially steady- 
state response characteristics of cascades of blades operating at this level 
of reduced frequency. Also, it was assumed that no crossflows or mixing 
occurs between sectors, which according to Reid (Reference 7) may be a good 
assumption for circumferential distortion patterns and according to Plourde 
and Stennlng (Reference 8) may be a good assumption for compressors with low 
gap-to-radius ratios. Finally provisions were made to determine the circum- 
ferential displacement of a streamtube for each sector of distorted flow 
through the compressor and the displacement of the streamtube was assumed to 
be equal to the circumferential displacement of the sector. 

Compressor maps and the associated inter-stage flow characteristics for 
distorted inlet flow conditions representative of the distortion tests con- 
ducted on the two J85-13 engines at the Lewis Research Center were then 
generated. The total number of corrected speed lines generated with inlet 
distortion was twenty five and the circumferential extent of the distortions 
considered was limited to 180”. The types of inlet distortion simulated 
corresponded to the total-pressure distortion imposed on the Moss engine, 
and the total-temperature distortion and a combination of the two with the 
distorted temperature region circumferentially opposed to, coincident with, 
and 90® overlapped with the distorted pressure region lnq>osed on the Mehalic 

engine* 


The computer model predictions were once again verified based upon the 
comparisons achieved between the predicted maps and the corresponding experi- 
mental maps noted in References 4 and 5. Specifically, verification was 
based upon the accuracy of predicting the flow-pressure ratio characteristics 
for each speed, and the surge lines. Following this verification, the pre- 
dicted stage velocity diagrams, blade and stator incidence angles, diffusion 
factors, and loss coefficients at surge and at two points below surge for 
corrected speed were analytically determined. 


Included in this report are: (1) The loss in surge pressure ratio at 

constant corrected speed and at constant corrected flow; (2) The ^plitude 
of total-pressure, static-pressure, and total-temperature distortion at the 
inlet, at each stage, and at the compressor exit; (3) The rotation (circum- 
ferential displacement) of the distorted sectors across each stage; and 
(4) The velocity diagrams, pressures and temperatures, and blade and stator 
diffusion factors for each stage for at least two sectors, depending on the 
type of distortion and its circumferential profile. Only the tables and 
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figures which serve to substantiate a point or which suimarizes final results 
are Included in the text. All documentation of model input and other addi- 
tional information is relegated to the appendix. 

The results are presented in the International System of Units. The 
dimensions of the compressor system which was designed using U.S. customary 
units are presented in both systems of units. 
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3.0 BLADE ROW DYNAMIC MODEL 


In this section is described the manner in which the Djrnainlc Digital 
Blade Row Compression Conqponent Stability Model is applied to the compressor 
of the J85-13 turbojet engine, the determination of the blade row character- 
istics, and the solution of the governing differential equations using a 
time-marching technique based on a Taylor series expansion. In addition, a 
steady-state model option is described which is used for determining some 
model input parameters and for initialization of the dynamic program. Also, 
the details of the parallel compressor option of the program are discussed. 


3.1 GENERAL DESCRIPTION OF ENGINE AND MODEL 

3.1.1 J85-13 Compressor 

The compressor of a J85-13 engine has eight stages with a variable camber 
IGV and variable third, fourth, and fifth stage bleeds located in each stator 
channel at the casing wall. The IGV trailing edge flaps and the bleeds are 
ganged together and are scheduled as a function of corrected speed biased by 
compressor-face total temperature. The geometry of the J85-13 engine in the 
compressor and combustor regions is shown in Figure 1. The nomonal IGV and 
bleed schedules are given in Figure 2. 

It is ^propria te to discuss the two engines modeled in this study. 

The pure total-pressure distortion patterns that were simulated during 
study were obtained during testing of a J85-13 engine known as the "Moss" 
engine (Reference 4) . This engine was run in support of the NASA casing 
treatment program. The clean inlet and distortion data utilized in this 
model were obtained from the untreated configuration with solid compressor- 
case inserts. The pure total— temperature distortion patterns and the com- 
bined total-pressure and total-temperature distortion patterns that were 
simulated were obtained during testing of a J85-13 engine kno%m as the 
"Mehallc" engine (Reference 5). 

The "Moss" engine was instrumented at the engine face with a 60 probe 
array (5 rings/ 12 rakes) to measure the total-pressure distortion patterns 
generated during testing. It was from these data supplied by the NASA Lewis 
Research Center that the patterns used in conducting the parallel compressor 
efforts of this program were selected. Further details concerning the in- 
strumentation can be found in Reference 4. The "Mehallc" engine was instru- 
mented at the engine face with a 60 probe array (5 rings/12 rakes), for 
nmasurlng total- temperature distortion and with an S^probe array for moni- 
toring the gross levels of total-pressure distortion. Further details can be 
found in Reference 5. 
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Figure 1 • Section View of J85~ 13 Compressor • 







3»1»2 Compressor Model 


For the purposes of this study, the compressor model Includes volume*! 

distortion measurement plane and do:ltre^"f 
(A?> ® f f ® located at the turbine nozzle diaphragm 

r^ili T*"® including these extra volumes is to Insure that 

realistic boundary conditions can be imposed. 

compressor model consists of twenty-nine volumes. There are 18 
bladed voltes (one blade row per volume) consisting of the IGV rotors 1-8 
stators 1-8, and the OGV and 11 free volumes. ThesI free vol^;s coLLt !f 
^o volumes be^een the instrumentation plane and the leading edge of Jhe 

Jot^r 1 aT 1-sding edgfof 

k" between the trailing edge of the OGV and the tur- 

bine diaphragm. The axial lengths of the free volumes are chosen to be c^ 
mensurate with the length of the longest blade chord axial XectlL^ afft 
if control the upper bound of the frequency response 

SL if -"ends fZ the trailing 

upstream stator to the leading edge of the downstream stator 

row extendi f^ii inter-blade row gaps while the length of a stator blade 

statir leading edge of the stator to the trailing edge of the 


in the model as well as the boundary layer blockage and 

hit r i '""a "" « 1 * the 02* pari^ters aS t“ 

blade leading edge and trailing edge metal angles, respectively. The given 

iLT^TaZl exit air angles. The Alheren^e^elntZ 

exit air angle and the metal angle is the deviation angle which for stators 

e? “ 
countered by the J85-13 compressors in this study. 

s,>ec±nlt^tl "" accomplished in a distributed, but un- 

Jakf fiii: Z Tnl ^ «« assumed to 

p ce in rotors, that is, no losses are accounted for in blade free 

^d^^^ ^ stator blade volumes, although this is not a restriction of the 

The r'nt deviation angles vary as functions of the Incidence angle 

Mehalic engines as a function of IGV Stator Setting (SS). 



” Moss 

M 

II 

Mehalic" 


SS 

lAA 

SS 

lAA 

80 

33 . 0 ° 

16 . 73 ° 




87 

16.5 

9.27 

6.7° 

3.98 

94 

1.0 

.61 

0.0 

0.0 

100 

0.0 

0.0 

0.0 

0.0 
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Table 1. J83-13 Compressor Geometry. 
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The reason that the above schedules are not functions of Tt 2 as in Figure 2 
is that the Tx2 bias was subverted during NASA engine testing in a manner 
that controlled the IGV and blade schedules to the 294,3®K (70® F) or higher 
temperature curve. The *’Mehalic” engine had a modified schedule. The actual 
schedules used during this study are given in Figure 4. 

At this point, it is worth mentioning that because a blade row formu- 
lation is being used, stage characteristics information was used in the form 
of a relative total-pressure loss coefficient and deviation angle rather than 
the more often used non-dimensional work and pressure coefficient as a function 
of flow coefficient. This approach tends to de-couple the inlet and exit 
stations of a blade row in a dynamic analysis since volume storage of mass, 
momentum, and energy are permitted within the work producing volume. Further, 
splitting a stage volume into two blade row volumes will potentially double 
the frequency capabilities of a dynamic compression component model. 

The third, fourth, and fifth stage bleed flows in the model are removed 
at the exit of the stator while holding exit air angle constant. The percent- 
age of inlet flow that is removed from each stage is given in Table 2. 


Table 2, Bleed Removal Schedule (Percent 



of Inlet 

Physical 

Flow) . 


Engine 

Corrected 

Speed 

Stage 3 

Stage 4 

Stage ! 

Moss 

80 

3.95% 

4.90% 

5.78% 

Moss 

87 

2.28 

2.83 

3.33 

Moss 

94 

0 

0 

0 

Moss 

100 

0 

0 

0 

Mehalic 

87 

1.20 

1.49 

1.76 

Mehalic 

94 

0 

0 

0 

Mehalic 

100 

0 

0 

0 
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The boundary conditions imposed upon the model consists of specifying 
the total pressure and total temperature entering volume 1 and specifying 
the value of the exit flow fxmction (Wv^/Ps) at the exit of the last volume 
and its time rate of change. Throttling along all speed lines was accomplished 
by decreasing the flow function at a rate of 7.5 units per second. This rate 
was established by throttling at a number of rates and selecting one which 
would not cause the throttled speed line to deviate from the speed line 
obtained from steady-state solutions. This rate is unique to the J85-13 
model and would have to be determined for each new compression component that 
is modeled. 


3. 2 BLADE ROW CHARACTERISTICS DETERMINATION 


As in any inter-blade- row compressor analysis, the performance of each 
blade- row or stage is described by a set of relations known as characteristics, 
which describe the manner in which work is input and losses are generated as a 
function of inlet conditions. As part of this program, it was necessary to 
generate a set of characteristics for both the "Moss” and "Mehalic" J85-13 
engines. The data used for obtaining the characteristics resulted from NASA 
tests of these two engines and were supplied by NASA. Generation of the non- 
dimensional characteristcs (flow, work and pressure coefficient) was accom- 
plished by compressor design personnel familiar in detail with the aerodynamic 
design of the J85-13 compressor using the General Electric Stage Character- 
istics and Stage Stacking computer programs. These programs provided a 
number of tecniques for determining stage characteristics from test data and 
calculated overall performance from a stage stacking of such characteristics. 
The techniques employed in this program for characteristics determination are 
discussed in the following paragraphs. Assumptions made in the analyses, as 
noted previously, include associating all the stage losses with the rotor and 
specifying the stator deviation angle to be constant over the whole compressor 
map. 


In the case of the "Moss" engine, instrumentation provided total- 
temperature and total- pres sure test data at the leading edge plane of each 
stator blade row at three radial immersions. Based upon radially area 
averaged test data (20% hub, 60% midspan, and 20% tip), initial non-dimensional 
stage characteristics calculations indicated certain inconsistencies in the 
data* These inconsistencies appeared in the form of calculated negative loss 
coefficients and efficiencies greater than one. In order to determine the 
possible sources of error, radial and axial distributions of total pressure 
and total temperature were plotted for each speed for at least the lowest and 
highest operating pressure ratios. Wall static pressures were also plotted 
as a check on the tip total-pressure levels. 

Two obvious problem areas were revealed by these plots. At 100 percent 
corrected speed, the stage 2 radial profile of total pressure was inconsistent 
with the profiles of adjacent stages and other speeds; that is, the pitchline 
value was considerably lower than hub and tip values. The tip and hub values 
were held constant and a curve was forced through them. Based on the shapes 
of the stages 1 and 3 radial profiles, a new mass-weighted average pressure was 
calculated for rotor 2 only as all other pressures were left unchanged. 
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The temperature measurements appeared inconsistent in several stages* 
particularly 2, 6» and 7« These data presented a dilemma since only two 
options were available in the current data reduction program - one uses 
measured temperatures at every stage and the other which uses no measured 
Internal temperatures. Since the discharge pressures and temperatures were 
consistent with data from other J83 tests, the latter option was used. 

This option calculates an overall polytropic exponent based on compressor 
inlet - and discharge - total pressures and total temperatures according to 
the equation: 


EX = In (Tj^/T-p^)/ ^ (P‘P 3 ^^T 2 ^ 


Ol) 


where 


Y 


2 ^ 

^Poly 


and npoly Is the poly tropic efficiency. This exponent was then used in 
conjunction with the average interstage total pressures to calculate the 
Interstage total temperature according to the equation: 



02 ) 


The advantage of maintaining a constant exponent rather than a constant 
poly tropic efficiency is that an iteration on y * f(T) is avoided. However, 
a variable y is used in the remaining calculations. 

Upon performing a cons tan t~speed stage-stack analysis it was observed 
that the non-dimensional stage characteristics did not adequately reproduce 
the speed line indicated by the test data. Consequently, adjustments were 
made to the characteristics in order to correct obvious stage mismatches and 
in the process to obtain a better match with the test data. Figures 5 and 6 
indicate the data points and the final non-dimensional stage characteristics 
used for the 100% corrected speed line of the "Moss'* engine. Stage 2 data 
were modified as previously discussed and this forced changes in the Stage 3 
characteristics. These changes represent the worst case in terms of amount of 
adjustment needed. The complete sets of non-dimensional stage characteristics 
for the "Moss" engine are documented in Appendix A. 

The in^lied relative total- pres sure loss coefficients and deviation angles 
were calculated along each speed line from the final set of nondimens ional 
characteristics. The relative total-pressure loss coefficient and deviation 
angle distributions for each rotor on each speed line were then curve fit with 
a least-squares polynominal in order to provide explicit analytical expressions 
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100% N/^^ . 


19 


PRESSURE COEFFICIENT 


1.20 


1.00 r 


■a- 


.80 r 


.60 


.40 


.20 I- 




LEGEND 


O ROTOR 

§ ROTOR 
ROTOR 
ROTOR 
ROTOR 
Di ROTOR 
Q ROTOR 
O ROTOR 


1 

2 

3 

4 

5 

6 

7 

8 


o I — — 1. — I » 

.40 .50 .60 .70 .80 

FLOW COEFFICIENT, 


Figure 6. Preliminary and Final Pressure Coefficients, "Moss" Engine 
10094 N/y/0. 


20 


Elu'. . : 

ora^L... 


I 


and avoid inefficient interpolation schemes in the dynamic model. Each 
characteristic was represented by a second, third or fourth order polynomial, 
whichever gave the best fit of the data points. In general, characteristics 
between the last data point and to some arbitrary point beyond the surge line 
were provided by extrapolation of the polynomial beyond which linear extra- 
polation was used as illustrated in Figure 7. Figures 8 and 9 illustrate the 
resultant distributions of relative total-pressure loss coefficient and de- 
viation angle used for the **Moss” engine 100% corrected speed line. Curves for 
the remaining speeds are given in Appendix A. 

The nondimens ional characteristics for the "Mehalic” engine were derived 
in a manner similar to that previously discussed. However, the "Mehalic** 
engine instrumentation provided casing static pressures as the only source of 
interstage data. This constraint necessitated the use of the constant poly- 
tropic exponent method of calculation in conjunction with casing static 
pressures in order to determine the stage performance and characteristics. 

These characteristics are documented in Appendix A. Further, Appendix A 
provides a compilation of the relative total-pressure- loss coefficient and 
deviation- angle distributions used for all speeds Investigated plus tabular 
listings of the coefficients of the polynomial representations. 


3. 3 ANALYTICAL TECHNIQUE 

3.3.1 Equations of Change 

The complete set of non-linear partial differential equations which 
describe the transfer and storage of mass, momentum, and energy within a 
fluid are called the equations of change (Reference 9). These equations 
have been Integrated once over an arbitrary volume of the flow system to 
obtain the macroscopic balances for quasi one- dimens ional flow without heat 
transfer and are reproduced below in the form in which they are used in the 
dynamic compression con^onent model. 
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The subscripted variables on the right-hand side of the equations refer to 
quantities at the inlet and exit of the control volume. Variables on the 
left-hand side refer to volume averaged quantities, l.e., in generalized 
form 


7 . J g dV 
^ ' / dV 
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The energy equation (Equation 3-5) was derived by combining the equation 
of change for energy and one of the thermodynamic TdS relationships. 


3.3.2 Force, Pressure and Entropy Production Terms 

This set of equations (other than being applicable to quasi one- 
dimensional flows without heat transfer and a finite, but small volume) 
properly and exactly describes the state of a fluid in motion. In order 
for a solution to be obtained, it is necessary to supply the caloric and 
thermal equations of state and expressions for Fg, and Sp*. 

Fb (Equation 3-4) represents the blade force acting upon the fluid. The 
blade force can be determined through reference to the following sketch: 
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and Is derived from the Euler Turbine Equation. The direction of the lift 
vector Is assumed to be 


= 2 <^i + ^2^ 


(3-9) 


The drag force (Fp) is obtained from the following equation 
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(3-10) 


which Is based upon an analogy with the drag coefficient for duct flows. The 
term Fd 2 in Equation 3-7 is then obtained from the relation Fj)^ “ Fjj/cos Boo. 


The prime (^) symbol Indicates the value of the parameters with respect 
to the relative velocity frame of reference and is the flow area per- 
pendicular to the direction of the entrance relative velocity vector. It 
should be noted that in steady flow a momentum balance, in general, will 
not give the same total-pressure rise per stage as does an energy balance. 
The reason for this difference is that the direction of the blade lift 
vector is not exactly the arithmetic average (Equation 3-9) of the flow 
angles. Comparison of the steady-state momentum and energy balance 
solutions permits the determination of a small "correction angle" which can 
then be added to to give the proper lift direction. It should be noted 
that the more familiar expression for the lift direction given by 


tan B« = 1/2 (tan + tan B 2 ) (3-11) 

also required use of a correction angle. Hence, Equation 3-9 is used in 
this formulation since it results in simple analytical relationships without 
introducing any compromises to accuracy or frequency response. Appendix B 
provides documentation of the correction angle relations used in this program. 
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(Equation 3-4) represents the mean pressure over the lateral surface 
area of the volume element as sketched below: 



where Pj^ = / PdA/ /dA. Although an analytical expression for the mean 
pressure acting on lossless, blade free volumes in steady flow can be 
derived, it leads to redundancy in the system of equations describing un- 
steady flow. Therefore, based upon steady-state momentum-balance analyses, 
an approximate linear expression for calculating blade-free volume mean 
pressure as a function of area convergence and inlet and exit pressures has 
been established for the J85-13 compressor model. Figure 10 illustrates the 
correlation developed in terms of a scale factor (FV) for zero-swlrl free 
volumes; its form is given by the following equation: 


Pm = ^1 + FV (P^) / (1 + FV) 


(3-12) 


Similarly, an additional correlation was established for the non-zero 
swirl free volume between the IGV and the first rotor. As shown in Figure 11, 
this correlation is a function of IGV exit air angle. It should be noted 
that the correlations presented in Figures 10 and 11 are unique to this formu- 
lation and the J85-13 compressor; they are probably valid only for other 
compressors with similar Mach numbers and area changes. 

For blade-row volumes, investigations have revealed that a good approxi- 
mation for the mean pressure is two-thirds the higher of the inlet or exit 
static pressures plus one third the lower pressure. Deviations from this 
approximation are accounted for in the lift direction correction angle. 
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Figure 11. Mean Pressure Correlation For IGV Trailing Edge Volume 



Sp (Equation 3-5) is the term which represents the total rate of 
irreversible conversion of mechanical to internal energy and, in the case 
of this model, represents the entropy production due to blade row losses. 
It can be obtained from the expression; 


S„ = WR In 
Jt* 




ideal 

actual 


(3-13) 


where the ideal relative total- pres sure ratio which accounts for the change 
in pitch line radius (Reference 10) from the entrance of a rotor blade row 
to its exit is written as 


!l2\ 

^ ' ideal 



1 



.t/(y - 1) 


(3-14) 


Mx is equal to the ratio of the blade row exit pitch line wheel speed to the 
inlet relative stagnation velocity of sound (2irNr2/aj^) . In the case of 
a stator, the ideal relative total-pressure ratio is equal to one. The actual 
relative total-pressure ratio requires knowledge of the relative total-pressure 
loss coefficient which is defined as 


**^2) " ^12) 

^ ideal ^actual 

la - 

PTi - Pi 


(3-15) 


Equation 3-15 can be rewritten in the form 



actual 



'ideal 


- 0 ) 



y/(y - 1) 


(3-16) 


Hence, Equations 3-14 and 3-16 when substituted into Equation 3-13 provide 
complete definition of the entropy production term. 


30 



As might be expected, the input to the program, in addition to physical 
speed, inlet conditions and conQ)ressor geometry, requires the relative 
total-pressure loss coefficient (to^) for the rotors, the tangent of the 
deviation angle (6) for the rotors, and the tangent of the correction angle 
(3p) for both the rotors and stators. The relative total-pressure loss coef- 
ficient and the deviation angle are derived from stage stacking results 
based upon clean inlet flow test data obtained by throttling at constant 
speed. The correction angle is obtained by comparing steady flow force and 
energy balance solutions on a blade row basis. These parameters can be 
represented as functions of incidence angle and are input to the program in 
this manner. These parameters, in conjunction with the velocity triangles 
and other ancillary relations, permit the determination of the thermodynamics 
of the fluid at each station. Blade-free volumes are treated as lossless 
volumes with no imposed blade force; hence, the Fg and Sp terms of Equations 
3-4 and 3-5 are identically zero. 


3.3.3 Calculation Technique 

Time dependent solution of the system of equations (Equations 3-3 
through 3-5 and the relations for Fg, Pj^, and Sp) that comprise the dynamic 
digital compression component model is effected through a Taylor series which 
establishes the values of the three independent volume- averaged variables at 
the next increment in time. In the case of this model and with reference to 
th e left hand side of Eqiiations 3-3, 3-4, and 3-5, the variables p, W, and 
pS are the ones for which a solution Is sought. Solution is now straight- 
forward and will be illustrated for one variable - the volume-averaged density. 
Considering that this method is applicable to any volume, the subscript 
indicating the k-th volume will be dropped. The Taylor series for volume- 
averaged density correct to second order can be written as: 

5(t+4t) - 5(t) + ^ 4t + ^ ( 3 - 17 ) 

3 1 


whore: 

p(t) is established by the initial conditions or from the previous time step. 


“ y (w^ ~ Equation 3-3 and differentiating 


Equation 3^3 with respect to time yields; 


gj.2 V \ at at / 


at 


(3-18) 



Examination of Equation 3-18 reveals that the right-hand side is composed 
of derivatives of station values of flow with respect to time. Since Equation 
3-4 will supply only the derivatives of the volume-averaged flow with respect 
to time, use of an interpolation scheme for obtaining station values from 
Volume— averaged values will permit Equation 3—18 to be solved for the second 
partial derivative of volume- averaged density with respect to time. Equations 
3-3 and 3-18 then can be substituted into Equation 3-17 to obtain the estimate 
of the volume averaged density correct to second order at the next Increment 
in time. Equation 3-18 implies that first derivatives with respect to time of 
a large number of terms (e.g. Fg, P| 4 , and Sp) will be required. Although 
these expansions are lengthy, they can be derived in a straightforward manner 
and will not be reproduced here due to lack of space. Similarly, this technique 
can be used for the remaining two variables (W and pS) and can be continued 
from one time step to the next for the desired number of time steps. 

This calculational scheme is numerically stable, and (to date) anomalous 
behavior has occurred only when a physical aerodynamic instability in the 
flow would be expected to occur. 

Now the calculational technique utilized in the Dynamic Digital Blade 
Row Compression Component Stability Model can be discussed and is illustrated 
in Figure 12 in block diagram format. The use of the blade row building 
block concept allowed construction of a generalized model which is indepen- 
dent of the particular compression component being simulated. 

Block I is a statement of the required dependent variable information, 
that is, volume-averaged density, flow, and entropy which are available from 
either a steady-state (SS) initialization or a previous time step of a time 
dependent (TD) analysis. 

Block II presents the macrobalances in the form they are used in the 
analysis. The variables on the right-hand- side of the equations are 
station-value properties. Knowledge of these parameters allows the first 
time derivatives of the volume averaged properties to be calculated. However, 
as stated in Block I, only volume- averaged quantities are available at the 
beginning of each time step. Therefore, Block II illustrates that it is 
necessary to interpolate between volume-averaged parameters in order to 
obtain station- value properties. 

In the case of blade-free volumes, where no blade forces or entropy 
production takes place, it is only necessary to calculate station axial 
velocity in order to evaluate the equations of change. As shown in the lower 
branch of Block IV, the assumption of constant absolute flow angle across 
the volume is made. Total pressure, total temperature and other desirable 
parameters are also calculated at this point. A special case of the blade 
free-volume calculations is the imposition of the boundary conditions. At the 
model Inlet, constant total pressure and total temperature as well as constant 
entropy are maintained. At the exit of a model of a turbojet compressor and 
burner, a specified exit flow function boundary condition is imposed. This 
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boundary condition was derived from the assumption that a fictitious, zero- 
length choked nozzle existed at the burner exit. The fictitious, zero- 
length choked exit condition implies the upstream flow function at the 
burner exit is specified and can be either held constant or changed at some 
rate to simulate a constant-speed throttling process. 

As indicated in the upper branch of Block IV, the presence of a bladed 
volume requires the net axial blade force, entropy production, and station 
axial velocities to be calculated. Calculation of the net axial force and 
entropy production terms require knowledge of the loss coefficient and devi- 
ation angle. This information is available as polynomial representations 
which are functions of incidence angle. Stationary blade rows are assumed to 
be lossless with constant deviation angles. 

Once the flow conditions at the stations are completely described, various 
quantities of interest can be calculated such as stage coefficients, diffu- 
sion factors, etc. With all the necessary quantities on the right hand side 
of the macrobalances available, the first time derivatives of the volume 
averaged properties can be calculated as indicated in Block V using the macro- 
balances of Block II. 

Expressions for the second time derivatives of the volume averaged 
quantities are obtained from differentiating the macrobalances with respect 
to time. Analytical expressions for the time derivatives of the station 
properties can be evaluated by interpolating between volumes and through use 
of the macrobalances. Thus, as indicated in Block VI, the second time 
derivatives of volume averaged properties can be calculated. 

This procedure for calculating station properties and evaluating first 
and second time derivatives of the volume averaged properties can be carried 
out for any number and types of volumes (Block VII) and is not dependent on 
the particular geometry being modeled. Once these calculations are carried 
out for all the volumes, the solution can be advanced to the next time step 
through use of the second-order Taylor’s series approximations. Block VIII- 
As specified in Block IX, the technique can be repeated for as many time 
steps as required by the event being simulated. 


3 . 4 STEADY-STATE OPTION - INITIALIZATION 


Although the dynamic compression system model is of prime interest In 
this study, many of the preliminary tasks required the use of steady-state 
analyses. In addition, all explicit time- dependent calculation schemes 
require Initial conditions with which to start the solution procedure. There- 
fore, an optional steady-state solution technique was included as an integral 
part of the computer program which can be used to supply initial conditions to 
a transient solution or to perform separate steady-state analyses. The steady- 
state option was made compatible with the equations of change of Section 3.3 
in that the time derivatives are set to zero and the downstream conditions at 
each volume are solved by satisfying the steady-state continuity and momentum 
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equations rather than the more often used energy equation. Con^atibility 
between the steady*~state and dynamic portions of the program was further 
maintained through the use of identical techniques for Including the loss 
coefficients, deviation angles, lift direction correction angles C&-). and 
the mean pressures (Pjj) . 


3.5 PARALLEL COMPRESSOR 


As it was the primary purpose of this program to analyze circumferential 
distortion, the macrobalances and the calculation technique were formulated 
in a manner such Chat a dynamic parallel-compressor analysis, constrained by 
the well-known parallel-conq)ressor boundary conditions, could be performed. 
That is, it was assumed that the conq>ressor was made up of several sectors, 
acting as Independent compressors operating with different inlet total- 
pressure and/or total- temperature conditions. Each sector operated on its 
clean inlet characteristics and exited to a common static pressure downstream 
of the compressor discharge diffuser. This exit point provided the only 
location where Inter^communicatlon between sectors was allowed; downstream 
of this point the flow was assumed to be completely mixed and was modeled 
as a single duct flow which terminated in a choke plane at the turbine 
diaphragm (A 4 ). The model developed during this study was unique in that 
the overall physical inlet flow was specified and the flow split to the 
various sectors was determined as a function of the imposed distortion levels 
and the calculated compressor diffuser-discharge uniform static-pressure 
boundary condition. 

When simulating a temperature distortion, the model incorporated an 
interpolation technique to determine blade-row characteristics for corrected 
speeds for which data were not available. The technique is illustrated in 
Figure 13. For each blade row at each corrected speed for which data were 
supplied, values of the characteristic at Intervals of constant percent 
tan(i) of the data range were established. The sets of characteristic data 
at constant percent intervals were then curve fit with a cubic spline to 
establish the blade-row characteristic as a function of corrected speed. 
Interpolating the spline— fit data at each constant percent interval for a 
particular corrected speed supplied a set of characteristic data as a function 
of tan(i). These data were then fit with a least— squares polynomial to pro- 
vide an analytical expression for the blade-row characteristic. Beyond the 
range of available data, the spline-fit curve fits were linearly extrapolated. 
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4.0 RESULTS 


Apart from determining suitable blade row characteristics for both of the 
J85-13 engines considered in this program, the objective of this study was to 
establish the validity and capabilities of the computer model. This involved 
verification of the ability of the pitch-line analysis technique to accurately 
represent the performance of the compressor and the demonstration of the 
ability to predict the stability limit of the compressor for clean inlet flows. 
Once the capability of the model was established, the stability limit of the 
compressor was predicted when subjected to circumferential, 180® 1/rev total- 
pressure, total-temperature, and combined total-pressure and total-temperature 
distortions. Combined inlet distortion conditions included configurations 
where the total-pressure distortion was opposed, coincident, and 90® over- 
lapped with the total-temperature distortion. 


4.1 CLEAN INLET ANALYSES 


Clean inlet compressor maps for both versions of the J85-13 considered 
are presented in Figures 14 and 15. It should be noted that the "Moss" engine 
clean inlet data (Reference 4) used in this study were taken from NASA data 
readings 521-543 which were obtained after the distortion tests had been run. 
The "Mehalic" engine clean inlet data (Reference 5) were taken from NASA data 
readings 13-38. As can be seen, there is excellent agreement between the 
NASA test data and the speed lines generated from the previously determined 
blade-row characteristics. Included on the figures are the stability limits 
obtained for each speed line by dynamic throttling simulations. The throttling 
simulations were accomplished by specifying the decrease in the flow function 
at the exit of the model as a linear function of time. The choked exit 
boundary condition was never explicity calculated, but rather was handled by 
assuming a zero length choked exit nozzle existed. As such, a change in the 
exit area of a choke plane would impress a change in upstream flow function. 

The rate of change of the flow function with time was chosen low enough such 
that the dynamic solution did not deviate from the steady-state speed line. 

It was the intent of this portion of the study to illustrate that the stability 
limit of a compressor could be determined from the dynamic response of a quasi- 
steady-state representation of the speed line. Figure 16 provides a more 
detailed view of the throttling process for the 94% N//© '^Moss" engine speed 
line. It should be noticed that up to time step 5000 the decrease in flow is 
well behaved, that is, the flow change per each 1000 time-step increment is 
roughly the same. However, in the region of the experimentally determined 
surge line, the flow decreased about twice as much in 700 time steps as it 
did in the previous 1000 time steps even though the throttling rate was main- 
tained at the same value as used throughout the throttling process. This 
behavior is typical once the stable operating region is exceeded. Special 
note should be taken of time step 5400 as further discussion will be undertaken 
later in this section. 
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Figure l6. Detailed Throttling Representation of "Moss" Engine at 94% , 



In an effort to establish stability criteria, many of the traditional 
compressor performance parameters were investigated for indications of 
anomalous behavior during the throttling simulation, particularly in the 
experimentally determined surge region. Figures 17 to 26 present the time 
history plots of several variables for the "Mehalic" engine, lOOX corrected 
speed line near the surge point. These plots show the time histories of: 

Overall Compressor Total-Pressure and Total-Temperature Ratio, 

Tangent of Incidence Angle at Entrance to Rotors, 

Rotor Diffusion Factors, 

Stage Total-Pressure Ratio, 

Rotor Axial Velocity Ratio, 

Stage Total-Temperature Ratio, 

Rotor Flow Coefficients, 

Rotor Work Coefficients, and 

Rotor Pressure Coefficients, respectively. 

As can be seen, none of these parameters exhibit anomalous behavior as 
the solution progresses through the region of the experimentally determined 
surge point, that is, the time when the model pressure ratio and corrected 
flow equal the values at the experimentally determined surge line. As a 
result, the performance of several of the dynamic variables of the throttling 
simulation was investigated. Figure 27 presents the time history response 
of the ratio of the rotor-volume-flow time derivative to exit-volume-flow 
time derivative. As the compressor is throttled, the parameter exhibits 
either a constant value of one or a monotonlcly increasing value in the region 
of the experimentally determined instability. Since the level of the exit 
flow derivative is constrained by the imposed exit boimdary condition, the 
behavior of the flow derivative ratio Indicates that internal perturbations 
are amplified significantly as the model nears the surge point. It was 
observed that once large amplifications had been encountered, the failure of 
the dynamic solution was imminent as It would predict impossible values of 
flow variables such as density being less than zero. In an attempt to deter- 
mine the level of flow derivative ratio associated with instability, several 
test cases were run in which the compressor model was throttled to various 
levels of flow derivative ratio. The throttling process was then terminated 
in order to let the solution stabilize. It was discovered that once the flow 
derivative ratio had reached a level greater than 2 in all the rotor volumes, 
the solution was unstable and termination of the throttling process would not 
prevent the solution from progressing into the post-surge region and failing. 
Therefore, this stability criterion was adopted and used throughout the pro- 
gram. It was noted that, in all cases, this criterion resulted in surge 
occurring in the region where the speed line slope approached zero. 

Discussion is redirected to Figure 16. The dynamic solution data point 
indicated by time step 5400 is the point where the flow derivative ratio 
attained the value of two or greater, in all blade rows, and is also the point 
where the slope of the speed line goes through zero. 
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Figure 17. Stability Criterion Investigation - Overall Compressor Total -Pressure 
Ratio Versus Time, 
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Figure l8. Stability Criterion Investigation - Overall Compressor Total -Temperature 
Ratio Versus Time, 
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Figure 19. Stability Criterion Investigation - Tangent of Incidence Angle at Entrance 
to Rotors Versus Time , 
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Figure 21. Stability Criterion Investigation - Stage Total-Pressure Ratio Versus Time 
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Figure 23- Stability Criterion Investigation - Stage Total -Temperature Ratio Versus Time 
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Figure 24. Stability Criterion Investigation - Rotor Flow Coefficients Versus Time 
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Figure 25. Stability Criterion Investigation - Rotor Work Coefficients Versus Time 
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Figure 26 . Stability Criterion Investigation - Rotor Pressure Coefficients Versus Time 
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Figure 27* Stability Criterion Investigation - Stability Amplification Flinction Versus Time 


At this point, it was appropriate to question whether this stability 
criterion represented a limitation of the model or did it in fact represent 
the actual behavior of the J85-13 compressor in the region of the clean inlet 
surge line. Discussion of the speed-line zero-slope phenomena at surge with 
compressor designers Indicated that they had not seen any data which would 
dispute the fact, although in general, data sufficiently close to the surge 
line to define the slope going to zero are not usually obtained, especially, 
at high speeds (steep speed lines) . However, one instance of high speed data 
was recalled where much data was obtained and a nuid>er of surges were incurred 
while investigating the near surge behavior of a near vertical high-speed, 
speed line. The data showed that the slope of the speed line did go toward 
zero in a small region near the surge line with small changes (decreases) in 
flow. 


It has been pointed out that the slope of the speed line must go to zero 
at the surge line in the presence of a choked exit boundary condition (choked 
turbine diaphragm) which is representative of compressor operation in the 
speed ranges studied in this program. Goethert et al, (Reference 11) have 
pj’oposed a stability criterion that explains thxs behavior. *Oie criterion 
was derived by application of the continuity equation to the stage volume 
and using a linearized representation of the stage characteristic in the 
region of interest. For a single stage followed by a choked nozzle, the 
requirement for stability can be expressed as 


3Wi^2 


3PT 


i+2 


3W, 


3 Pip 


i+2 


(4-1) 


The first term represents the characteristic of the nozzle in terms of 
the manner in which the exit flow rate responds to a change in stage exit 
total pressure. The second term represents the stage characteristic in terms 
of the manner in which the stage inlet flow responds to a change in compressor 
exit total pressure. Equation 4—1 can be written in an equivalent form as 
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(4-2) 


where the subscript indicates the stage entrance conditions, the subscript 

indicates the choked nozzle throat conditions, and mp is the slope of the 
stage characteristic at the point under consideration. The first term is 
clearly positive and bounded. As one progresses along the speed line from 
high flow to low flow, the sign and magiitude of the second term is governed 
by the local slope mi. The slope is initially negative, goes to zero, and 
then becomes positive. In the region of zero slope, this term ejdiibits a 
discontinuous behavior as it goes from indeterminately large negative values 
to indeterminately large positive values for small changes in flow on the 
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order of e. It is at this point that the model will ejdiibit the character- 
istics of incipient instability since the sign of Equation (4-2) changes from a 
positive to a negative value. Hence, it was for this reason that it is 
expected that the speed line should exhibit zero slope at the surge line con- 
sidering that a choked turbine diaphragm boundary condition has been imposed. 
Tlierefore, it was possible to draw the inference that the more sophisticated 
solution of the flow aerodynamics represented by the subject Dynamic Digital 
Blade Row Compression Component Stability Model supports the stability 
criterion based upon a simplified, linearized model of the flow proposed by 
Goethert et al. 

Figure 28 presents the ratio of exit-to inlet-flow for the rotor volumes 
as a function of time. As the compressor is throttled toward instability, the 
data illustrate that the rotors exhibit flow storage and flow evacuation, that 
is, more fluid enters rotors 1 and 2 than leaves and more fluid leaves rotors 
3-8 than enters near the region of instability. If it is assumed that surge 
or overall compressor instability is associated with flow blockage, then this 
parameter could be useful in identifying the stage where the surge event is 
initiated. 

Sample documentation of the clean— inlet compressor performance is con- 
tained in Table 25 of Appendix C. 


4.2 TOTAL-PRESSURE DISTORTION ANALYSES 


The dynamic parallel-compressor model was used to analyze the response of 
the "Moss" J85-13 engine to 180°, 1/rev circumferential total-pressure distor- 
tions. Three different distortion-screen porosities were tested on the engine 
«ind the effects of total— pressure distortion produced by the screens were 
simulated by the model. Table 3 provides a tabulation of the speed lines 
investigated and the distortion levels imposed for the NASA data readings 
recorded nearest to surge. The total-pressure distortions were modeled as 
two, 180° sectors with levels ( (max-min) /avg) which ranged from 1.9% 

to 13.7%. Maximum and minimum values were established by averaging all probe 
readings in the high- and low-pressure regions respectively, disregarding 
probe readings which were determined to be erroneous through examination of 
distortion profile plots. Sample plots of the radial and circumferential 
profiles of the distortion data analyzed are given in Appendix D (Figures 80 
through 82). 
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Figure 28. Stability Criterion Investigation - Ratio of Volume Exit Flow to Inlet Flow 


Table 3, 180®, 1/Rev Total- Pressure Distortion Cases 


%N/Ve 

Screen* 

Near Surge 
Reading 

APp/Pr** 

80 

4M 

490 

0.0187 

80 

9M 

126 

0.0475 

87 

4M 

485 

0.0281 

87 

7-1/2M 

379 

0.0554 

87 

m 

90 

0.0677 

94 

4M 

481 

0.0362 

94 

7-1/2M 

381 

0,0751 

94 

9M 

94 

0.0967 

100 

4M 

478 

0.0524 

100 

7-1/2M 

384 

0.1119 

100 

9M 

99 

0. 1367 


* Percent Op>en Area 
4M - 74.0 
7-1/2M - 57.8 
9M - 50.7 


** = (P - P 1/p 

T T T Max T Min‘^^*T Avg 


n. H 1 ^ I? present the performance of the parallel-compressor 

r?/L^"oM r ’“!!! vf for the three distortion screens - 4M, 

^ Table 3). Shown on the figures are the NASA test data, experi- 

mental distortion surge line, experimental clean-inlet surge line the 
parallel-compressor model speed-lines, and the stability limit predicted by 
the model. Regions on the figures where the distortion surge line exceeds the 
s„rge line were felt to be the result ot using clean Inletdata 
stained at the end of the engine test and as such represents the petforaance 
of a degraded ™gi„e. In all cases, the distortion was slnnlated as two, 180- 
sectors with the distortion levels deterined fro. the near-surge deta Md 
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Figure 31. Effect of l80®, l/Rev Total-Pressure Distortion (9M Screen) on Surge Line 
of *'Moss** Engine. 



held constant over the entire speed line. The stability limit of the overall 
compressor was established as that point where the previously determined 
stability criterion was exceeded by any of the parallel-compressor sectors. 

For all distortion levels, the model exhibited the ability to predict the 
experimental surge point quite well at 100% corrected speed. Conversely at 
the lower speeds, the model predicts a more conservative loss in surge pres- 
sure-ratio than demonstrated experimentally. This result was expected as 
even with small distortions, the parallel-compressor concept maintains the 
Integrity of each sector and prevents any mixing or redistribution of flow 
Further discussion of this result is reserved for Paragraph 4.5.1. The 
apparent inability of the model to match the test data at 94% corrected speed 
for the higher distortion levels was not resolved. Since it was possible to 
match the test data at other speeds, erroneous test instrumentation was not 
felt to be a probable cause. A significant factor might possibly be the 
proximity of the 94% corrected speed to the break-point on the IGV and bleed 
schedules For the ■’Moss" engine (Figure 4) and the possibility of the distorted 

flow compressor not reacting to the clean- inlet- flow IGV and bleed schedules 
as expected. 

Sample documentation of the total-pressure distortion parallel-compressor 
analyses is given in Table 27. The predicted operating points for the indi- 
vidual parallel compressor sectors on the compressor maps are shown in 
Figures 83 through 85 of Appendix D. 


4.3 TOTAL-TEMPERATURE DDISTORTION ANALYSES 

A total of three 180“, 1/rev, circumferential total-temperature distortion 

"Mehallc" J85-13 engine. Distortion levels of 
ATt/Tt ((max-mlnUavg) from 3.6% to 15.6% were investigated. Table 4 pro- 
vides a compilation of the distortion levels simulated, the NASA test reading 
for the points recorded near surge from which the data were obtained, and the 
nominal corrected speed of the speed lines simulated. As shown in Figure 86 
of Appendix E, the profile plot of the circumferential total-temperature 
distortion distribution does not exhibit a very "square" profile. In addi- 
tion, the radial profile plot indicates the existence of substantial radial 
distortion as well. The procedure for modeling the temperature distortion as 
two, 180 sectors was to radially and circumferentially average the probe 
readings from the center four rakes in the high- and low- temperature regions 
and the resultant temperatures were then assumed to extend over their respec- 
tive 180 sectors. Any differences in corrected speed and total-temperature 
distortion levels noted between the results presented In this report and the 
results presented in Reference 5 are the result of the above manner of averag- 
ing probes in the high and low tempe^ture regions to obtain the maximum and 
minimum values to calculate the AT^/f^ parameter. Preliminary 4-sector model- 
ing efforts which provided for sectors of Intermediate temperature between 
the high and low regions yielded near-identical overall-average temperatures 
as the two-sector model. Thus both simulations would Indicate the same 
corrected speed, and overall performance, but the 4-sector analysis would 
require twice the computational time by requiring the consideration of two 
^ non-critical sectors. Bleed and IGV schedules were specified as 

a function of the average inlet total temperature. 
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Table 4, 180°, 1/Rev Total-Teraperature 

Distortion Cases. 



Near Surge 


%NA/0 

Reading 

ATjAt’ 

86.8 

568 

0. 1557 

93.4 

154 

0 . 0466 

99.6 

138 

0.0363 


- (T^ Min^^T Avg 


Figure 32 illustrates the results of the total-temperature distortion, 
parallel-compressor analyses of the Mehalic engine. The compressor map illu- 
strates the NASA test data, the experimental surge line, and the throttling 
simulation stability limits. A detailed view of the performance of the 
individual parallel compressor sectors is provided in Figure 87 of Appendix 
E. The experimental distortion surge line is shown disconnected as this indi- 
cates two levels of distortion amplitude (Table A) . A comparison of the test 
data and the model-generated speed lines indicates the validity of inter- 
polating blade-row characteristics to determine speed lines not established 
by engine tests. The deviation of the model speed line at 86.8 percent cor- 
rected speed was due primarily to the large extrapolation necessary to provide 
characteristics data for the high total-temperature sector operating at 83.6% 
corrected speed; an extrapolation equivalent to 15% of the data range beyond 
the lowest speed data available. As in the case of total— pressure distortion, 
the total-temperature distortion results indicate the model predictions corre- 
late quite well with experimental results at high corrected speed, but are 
conservative at the lower corrected speeds. At 93.4% corrected speed, the 
parallel-compressor results indicate that the associated sector at 92.3% 
corrected speed (See Figure 87 of Appendix E) did not experience instability 
until beyond the experimental clean inlet surge line. This was a result of 
only using data at 87, 94 and 100% corrected speeds for the basis of inter- 
polation; the inclusion of characteristic data at additional speeds would 
undoubtedly improve the ability of the sector to better match the experimental 
clean inlet flow surge line. In all the temperature distortion cases investi- 
gated, the low corrected speed sector (high T»p) was limiting and experienced 
the instability. Sample documentation of the total-temperature distortion 
analyses is contained in Table 28 of Appendix E. 
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”Mehalic** Engine. 


REPRODUCIBILrrY 0? THE 
v^UUllN'AL PACE IS POOR 


4.4 COMBINED TOTAL-PRESSURE AND TOTAL-TEMPERATURE DISTORTION ANALYSES 


In addition to total-pressure and total- temperature distortion analyses 
alone, analyses of the *’Mehalic” J85-13 engine were performed for combined 
total-pressure and total-temperature distortions of 180° in extent. Three 
orientations of the distortions were considered: opposed, coincident, and 

90° overlapped. The '^opposed pattern*’ term refers to the low ?j sector being 
opposed to the high Tx sector, while the "coincident pattern" term refers to 
the low Px sector overlapping the high Tx sector. With the exception of the 
90° overlapped configuration, the temperature distributions imposed on the 
model were determined in the same manner as discussed in the previous section. 
Total-pressure instrumentation at the upstream distortion instrumentation 
plane consisted of only eight probes (four rakes with probes at two immersions) 
as indicated in Reference 5. A square wave circumferential total-pressure 
profile was assumed. Examination of the limited test data indicated no signi- 
ficant radial total-pressure gradients and as such, the distortion levels 
imposed on the model were obtained from averages of the probe readings in the 
high- and low-pressure regions. In the case of combined distortions any 
differences in corrected speed and total-temperature distortion levels noted 
between the results presented in this report and the results presented in 
Reference 5 are the result of averaging the total- temperature probe readings 
as for the case of the pure total-temperature distortion patterns (Paragraph 
*rho dl ffcrcmcs between the total-pressure distortion levels reported 
in Reference 5 and herein are due to the difference in the distortion param- 
eter (Plavg ~ PTmin)/^Tavg Reference 5 versus the distortion parameter 

(PTmax “ ^Tmin) /PTavg ^sed in this report. Further, differences will arise 
due to the averaging of the total-pressure probe readings as noted above. 
Figures 88 through 90 of Appendix F contain sample plots of the circumferen- 
tial and radial profiles of the distortion patterns considered. 


4.4.1 Opposed Orientation 

Table 5 presents a listing of the opposed combined-distortion analyses 
performed, the distortion levels imposed, and the corresponding NASA data 
reading numbers. Since the distortions were opposed, it was possible to use 
a two, 180° sector parallel compressor model. Figures 33 and 34 present the 
results of the parallel-compressor modeling and the predicted stability 
limits for the cases of low- and moderate-temperature distortions, respec- 
tively. The stability-limit throttling simulations correctly predicted the 
qualitative change in the surge line with the different imposed distortion 
patterns. In addition, a comparison between the low-speed simulations of both 
figures indicated the model was able to predict the drop in the surge line as 
the level of temperature distortion was increased. As in the previous distor- 
tion cases, the model prediction of loss in surge pressure ratio corresponds 
quite well with experimental results at high speed but tends to over-predict 
losses at the lower speeds where the compressor speed lines are of lower 
slope. A comparison of the approximately 87% corrected-speed lines of Figures 
33 and 34 indicates the effect of extrapolating blade-row characteristics. 
Figure 33 demonstrates the model does a very credible job of matching the test 
data with 4.9% temperature distortion while Figure 34 illustrates a greater 
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CORRECTED FLOW (kg /sec) 

Figure 34. Effect of l80°, l/Rev Combined Total-Pressure and Total-Temperature Distortion, 
Opposed Orientation Moderate Total-Temperature Distortion on Surge Line of 
"Mehalic" Engine. 


deviation between the test data and model simulation with 9% total-temperature 
distortion. In both cases (87% corrected speed) the total-pressure distortion 
levels were equivalent. 


Table 5. 180®, 1/Rev Combined Total- Pressure 

and Total-Temperature Distortion Cases 
Opposed Orientation. 


%N/Ve 

Near Surge 
Reading 

87.2 

494 

92.5 

479 

86.8 

441 

99.4 

416 

* - <T^ „„ 

** Max 


ATiji /Tx * 

AP,p/Pj.*» 

0.0900 

0.0648 

0.0919 

0.0870 

0,0489 

0.0612 

0.0362 

0. 1262 


T )/T 
T Min T Avg 


p ) /p 

T Min^^^T Avg 


The limiting sector, in all the cases but one, was the low corrected- 
speed sector; the total-pressure distortion for the 100% corrected speed 
event was sufficiently high (12.6%) to cause the high corrected— speed sector 
to experience the limiting instability. Appendix F contains detailed compres 
sor maps Illustrating the performance of the individual parallel compressor 
sectors (Figures 91 and 92) . A comparison of the detailed maps for pure 
total-temperature distortion and the opposed combined distortions indicates 
the effect of total-pressure distortion in changing the operating points of 
tlie sectors. 


^.^.2 Coincident Orientation 

The second combined-distortion configuration considered was that where 
the total-pressure distortion (low Pf) and total- temperature (high Tx) are 
completely coincident (overlapped). The speed lines investigated, the dis- 
tortions imposed, and the corresponding NASA readings are listed in Table 6. 


67 



Appendix F contains sample plots of the circumferential- and radial-profiles 
for one distortion case as obtained from the NASA data. Figure 35 illustrates 
the resultant compressor map including the test data, parallel-compressor 
speed lines and the throttling simulation predicted stability limits. The 
experimental stall line is shown disconnected and represents two different 
levels of temperature distortion. As can be seen, the 2-sector parallel- 
compressor modeling of the distortions correlated well with test data and 
predicted the increased drop in the surge line due to the overlapping of the 
distortions. The detailed compressor map (Figure 93) contained in Appendix F 
illustrates the fact that the drop in surge line can be attributed to the 
effect of the parallel-compressor uniform-exit-static-pressure boundary 
condition in forcing the distorted sector further up its speed line than if it 
were operating with one distortion only. A comparison of similar operating 
pQifits for the coincident and opposed distortion xllustrates that while both 
parallel compressor simulations exit to the same pressure, the pressure dis- 
tortion sector in the coincident configuration must operate further up the 
speed line in order to produce sufficient pressure rise as it is operating at 
a lower corrected speed due to the superimposed temperature distortion. Again, 
it is apparent from Figure 35 that the model can provide an accurate estimate 
of loss in surge pressure in situations where speed lines are nearly vertical. 
At lower corrected speeds, where the lower slope speed lines require larger 
changes in flow for changes in pressure than at the higher speeds, the parallel 
compressor model forces the critical sector closer to surge than if the speed 
lines were more vertical, thus producing more pessimistic predictions of the 
surge point. 


Table 6. 180°, 1/Rev Combined Total- Pressure 

and Total-Temperature Distortion Cases 
Coincident Orientation. 


%NA/e 

Near Surge 
Reading 


ABj./P^** 

87.3 

504 

0.0873 

0,0675 

93.1 

436 

0,0371 

0.0914 

99, 1 

423 

0.0402 

0. 1272 

* AT /T = fT 

T ^ /T 



T T T Max 

T Min'' T Avg 

+♦ Ap /p ^p 

P VP 



T T T Max 

T Min ' T Avg 
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Figure 35» Effect of l8o°, l/Rev Combined Total-Pressure and Total- Temperature Distortion 
Overlapped Orientation on Surge Line of **Mehalic” Engine, 


4. 4,3 90*" Overlapped Orientation 


The final combined-distortion configuration investigated consisted of 
the condition where 180° extent total-temperature and total-pressure dis- 
tortions were overlapped by 90°. As presented in Table 7, total-temperature 
distortions ranged from 3.5 to 7.8 percent while total-pressure distortions 
ranged from 6.5 to 12.7 percent. The compressor response was analyzed using 
a four 90° sector, parallel-compressor model. Whereas total-pressure distor- 
tion levels were obtained as mentioned previously, total- temperature levels 
in each sector were obtained by averaging all inlet temperature probe data 
within the sector. Figures 36 and 37 present the results of the dynamic 
parallel-compressor modeling for low and moderate temperature distortions, 
respectively. Although the throttling simulations produce an accurate repre- 
sentation of the test data points and the qualitative changes in the surge 
line are duplicated, there is a greater deviation from the experimental surge 
line than with coincident distortions. This was due to the fact that the 
four-sector, 90° overlapped model (Figure 94 of Appendix F) had a higher exit 
static pressure than the coincident model (Figure 93 of Appendix F) and thus 

forced the critical sector further up the speed line producing a more pessi- 
mistic estimate of loss in surge pressure ratio. It is apparent that in 

order to make a reasonable analytical estimate of loss in surge pressure for 
small sector distortion, a more realistic model of the flow processes involved 
will have to be constructed. 


Table 7. 180®, 1/Rev Combined Total- Pressure 

and Total-Temperature Distortion Cases 
90° Overlapped Orientation, 


Near Surge 


%N/Ve 

Reading 

86.8 

445 

93.0 

432 

99.0 

420 

87.2 

499 

♦ AT /T 

= (T 

T T 

T Max 

AP /"P 
^T^ T 

= (P 

T Max 



APt/Pt** 

0.0412 

0.0651 

0.0346 

0.0888 

0.0352 

0.1272 

0.0781 

0.0672 


T )/T 

T Min ^ T Avg 


Avg 
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of ''Mehalic'* Engine. 


Sample documentation for one 90° overlapped case is given by Table 29 of 
Appendix F. Also contained in Appendix F is the detailed compressor map 
(Figure 95) for moderate total-temperature distortion. 


A. 5 SUMMARY OF PARALLEL COMPRESSOR RESULTS 

4.5.1 Distortion Sensitivity 

One of the most common techniques for indicating the response of a com- 
pressor or engine to various inlet conditions is to compare the surge pressure 
ratio of the distorted inlet-flow compressor to the clean inlet-flow compres- 
sor. This loss in surge pressure ratio (APRS) can be calculated on a constant 
speed or constant flow basis depending on the desired application of the 
results. As an aid in evaluating the overall ability of the dynamic parallel- 
compressor model to predict the degradation of the surge line when subjected 
to various inlet conditions, a comparison between the experimental- and 
predicted— loss in surge pressure ratios has been made. Figures 38 and 39 
respectively present the loss in constant speed and constant flow surge pres- 
sure ratios predicted by the parallel-compressor model compared to the experi- 
mental losses in surge pressure ratio.. At high corrected speeds, the parallel 
compressor model does a very credible job of predicting the actual loss in 
surge pressure ratio for 180 degree, 1/rev distortions. The data points on 
are scattered furthest from the perfect agreement line are 
indicative of results at lower corrected speeds and of the 90° overlapped 
distortions. It is felt the inability of the model to match the experimental 
performance is due largely to the parallel-compressor limitations which pre- 
vent any flow redistribution or communication between the sectors. 


Due to the close proximity of the blade rows to one another, it is felt 
that the major redistribution events probably have to take part in the free 
volumes upstream of the inlet guide vane. This line of reasoning is supported 
by a number of studies. Plourde and Stennlng (Reference 8) concluded that 
for compressors with normal clearances, circumferential flow redistribution 
occurs ahead of the compressor, and very little crossflow occurs within the 
compressor. Spring (Reference 12) extended the analytical technique of 
Plourde and Stenning (Reference 8) to determine the effect of the compressor 
pumping characteristic on the flow field entering the compressor. He found 
that significant static-pressure distortions and tangential velocity compo- 
nents existed at the compressor face while none existed far upstream where a 
pure total-pressure distortion was being impressed. For a pure sinusoidal 
circumferential total-pressure distortion of ±5% far upstream of the com- 
pressor, xt was found that at the compressor face, ±4% static pressure dis- 
tortH.ii was created with a concomitant tangential velocity distribution of 
17/. ol tlu- nN.an axial vcl.x ity. However, the corrected flow distortion 
drop|)ed Irom approximately 135% to approximately ±5%. Adamezyk (Reference 13) 
also sliows a flattening of axial velocity profile as the compressor face is 
approached from upstream. 
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Figure 39- Comparison Between Actual and Predicted Loss In Surge Pressure Ratio at Constant 



Braithwaite, Graber, and Mehalic (Reference 14) have derived a simplified 
parallel compressor model which yielded substantially the same results in 
terms of APRS) Calculated versus APRS)neasured correlation. Hieir analysis 
included all of the pure total-temperature and combined total-temperature and 
total-pressure distortion data points reported upon in this report. Further, 
it is noted that the Dynamic Blade Row Parallel Compressor model results also 
showed that opposed orientation total-pressure and total-temperature distor- 
tions can offset each other, thereby producing insignificant losses in surge 
pressure ratio. This is in agreement with the experimental data and the 
results of the simplified parallel compressor model proposed by Braithwaite, 
et al • 


4.5.2 Transmission of Distortion 

It is of interest to compressor design personnel to determine the con- 
tributions of individual stages to the attenuation or amplification of dis- 
tortion. The normalized circumferential total-pressure, total-temperature, 
and static-pressure distortions at the exit of each rotor have been calculated. 
Tlie distortion values were calculated as the difference between the highest 
and lowest circumferential values of the parameter normalized by the area 
y 0 ighted average of the parameter. As an example of the method of presenta- 
tion, the results for inlet total-pressure distortion at high speeds (94 and 
100%*N//0) are shown in Figures 40 through 42. Additional distortion ampli- 
fication results are given in Figures 96 through 110 of Appendix G. 

Examination of Figure 40 shows that the total-pressure distortion is 
amplified in rotors 3 and 4 at 94% N//e while it is amplified slightly or is 
transmitted with no change by rotor 4 at 100% N//e. All other stages at both 
speeds attenuate the total-pressure distortion. figure 41 shows that with an 
inlet total-pressure distortion at 94 and 100% N//0, little temperature dis- 
tQj*tion is created or amplified until the pressure distortion reaches the 
rear stages (6, 7, and 8). The static-pressure distortion which accompanies 
the total— pressure distortion behaves in the same manner as the total— pressure 
distortion. In fact, this observation holds true for all the distortions 
simulated during this study. 

It Is Interesting to note the behavior of the "Mehalic” engine when sub- 
jected to 180°, 1/rev total-temperature distortion (see Figures 99 through 
101). Although there was no imposed total-pressure distortion, significant 
total-pressure distortion was generated as a result of the sectors operating 
at; different corrected speeds. At 99.6% corrected speed, a small total- 
pressure distortion was created in rotor 1 and attenuated thereafter in the 
compressor. At lower corrected speeds, substantial total— pressure distortion 
was generated in rotors 2 through 4 although rotor 2 significantly attenuates 
this distortion at 93.4% N//e. The magnitude of the imposed total-temperature 
distortion changes little throughout the compressor at the higher speeds with 
some slight attenuation taking place at the 86.8% corrected speed condition. 
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Figure 4 l. Predicted Total -Temperature Amplification For Inlet Total-Pressure Distortion at 
94% and 100% N/\/9. 


SPEED 



79 


Figure 42. 


An examination of the predicted amplifications for the combined total- 
pressure and total-temperature distortion, opposed orientation ^^^^ures 
throu^ 104) reveals for all cases, that the imposed total-pr^sure distortira 
is attenuated through-out the compressor. Comparison of the two cases at t e 
low corrected speed Indicated that the higher total- temperature-distortion 
aided in attenuating the total-pressure distortion levels. In all cases there 
is a slight attenuation of the imposed total- temperature distortion. 

The response of the compressor to combined total-pressure and total- 
temperature ^incident distortion (Figures 105 throu^ 107) is 
from the opposed combined distortion. The imposed total-pressure distortion 
reacts in ^manner similar to the "Moss" engine with pure total-pressure 
distortion. At the lower speeds, rotors 2 through 4 amplify 

while at 99.1 percent corrected speed, only rotor 4 causes amplification. In 
contrast to the opposed distortion pattern orientation, the coincident Pattern 
ner^otll-Lmperaiure distortion at the compressor exit was amplified somewhat. 

The distortion level response to the combined total-pressure ^d 
temperature distortion 90’ overlapped orientation (Figures 108 through 1 ) 

paralleled the coincident pattern distortion cases. 

To further understanding and to improve predicting the response of a 
compressor to distortion, it would be desirable to have ^ esti^te of t^ 
ci«umferential growth of a distorted sector as it passed throu^ the com- 
pressor. In order to accurately predict such a phenomenon it would be 
Lcessary to calculate the flowfield in detail between blade rows. Particu- 
larly in the region of the interface between the sectors. A task this 
magnitude falls outside of the scope of parallel-compressor representation , 
bursome estimates of the rotation of the sectors and their potential overlap 
as they travel through the compressor have been made. Based “P°" averag 
flow-JIgle technique, the parallel-compressor 

overlapping of sectors is probably not a powerful miring force for this 
particular compressor since the sectors tend to overlap by less^^ 5 
TOSt. This value is arrived at by comparing values of the paraMter ^0 
given in the sector performance tables of Appendices E and F. It is ^eco^ 
Jized however, that this analysis has not taken into account the locally hi^ 
val^s ^ Induced swirl which can exist at the edges of a distortion pattern 

at the compressor face. 


4.5.3 Diffusion Factor Analysis 

The blade row model employed in this study provides 
about the inlet and exit velocity diagrams of a blade row make it poss 
to calculate one of the parameters useful to the compressor designer, the 
diffusion factor. As an example of this type of ^alysis. rotor diffusion 
factors are plotted for the "Moss" engines with clean inlet flow for 80 and 
100 percent corrected speeds in Figures 43 and 44, respective y. ranee 

speed the diffusion factor is given for three flows - unthrottled, mid range, 
aL near surge. The largest difference between the two speeds occurs in roto 
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1 where at 100 percent corrected speed, rotor 1 has a considerably higher level 
of loading. Both figures show that as surge is approached, the level of dif~ 
fusion factor increases except in rotor 3 at 80 percent corrected speed, where 
it remains essentially constant and at 100 percent corrected speed where rotor 
4 shows a small decrease in the diffusion factor as the near surge point is 
/ipproached . 

Analysis of the critical sector for the 180°, 1/rev total-pressure dis- 
tortions shows that the value of the diffusion factor at surge is the same 
as that at surge for clean inlet flow. This result is expected since the 
parallel compressor sectors operate on their clean inlet characteristics. 
Investigation of the diffusion factor behavior when 180°, 1/rev total-tempera- 
ture distortions were imposed was not conducted since the critical sector 
operated at a different corrected speed than the average corrected speed of 
all sectors, the only speed for which the clean inlet surge line diffusion 
factor data was available. 


4.5.4 Critical Stage Analysis 

In conjunction with the development of a stability criterion discussed in 
Paragraph 4.1, the use of the parameter Wi+^/Wi (the ratio of the physical 
flow leaving a volume to the physical flow entering a volume) as an indicator 

of the blade row location of the first aerodynamic instability initiation was 
presen ted . 

Although it is felt that more experience concerning the behavior of this 
parameter must be gained before confidence in our ability to interpret its 
results is fully established, the parameter has been examined for the clean 
inlet flow and 180°, 1/rev inlet total-pressure distortion cases. The results 
of this analysis are given in Table 8. The indicated rotor is the one which 
shows Wi+i/Wj to drop first and also to give the lowest value of W^.j/Wi at 
the time the stability criterion has a value of 2 or greater in all the blade 
rows. Examination of these clean inlet flow results show that the blade row 
characteristics derived for the 'Hoss" J85-13 engine imply that flow breakdown 
IS initiated in the front of the compressor (stages 2 and 3) and does show a 
slight rearward movement with increasing corrected speed as would be expected. 
Since the blade rows in each sector of the parallel compressor are constrained 
to operate on blade characteristics determined from clean inlet flow data it 
is anticipated that the critical sector blade row where flow breakdown is’ 
initiated with inlet total-pressure distortion would be the same as the clean 
inlet flow results. The parallel compressor results also shown in Table 8 
substantiate this point of view. 
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Table 8. "Moss” 


Engine Stall Site Analysis. 


Percent 


180®, 

1/Rev Total- Pressure 

Distortion 

Corrected 

Speed 

Clean 
Inlet Flow 

Low 

Moderate 

Hi^h 

80 

2* 

2 

— 

2 

87 

2 

2 

2 

2 

94 

2 

2 

2 

2 

100 

3 

3 

3 

3 


* Stage where aerodynamic instability originates. 


Wenzel, Moss, and Mehalic (Reference 5) have determined the location of 
the "stall sites", from an analysis of the high-response interstage static- 
pressure data obtained during testing of the "Moss" engine. The results 
shovm below indicate that the "stall sites" are located in the rear half of 
the compressor for both clean inlet and distorted inlet flows. 

Stall Sites 

% N//e Clean 180®, 1/rev Pj 


80 5 ? 

87 - 5/6 

94 6/8 6 

100 7 6 

The model results and the interpretation of the test data do not indicate 
the type of agreement one would like to see. However, this lack of agreement 
may not be due to model deficiencies. Before judgment can be rendered on the 
ability of the model to predict stall sites, it is necessary that a compres- 
sion component which had extensive interstage instrumentation during testing 
be modeled and the latest techniques for analyzing and interpreting interstage 
data be applied. The problem of determining the location of hub stall sites. 
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‘l”ce It is unknown how long It takes for a 
hub stall pulse to ptop.gate to the wall easing «,d be detected bfuLl 
static-pressure transducers. ^ 

P^^ediction Of the blade row where flow breakdown is 

^haJ ri r" ^‘^‘=°”'Plished by a model, but at the same time it is clear 
that this is an area which requires a lengthy study in itself if oronar 
intetpretatlon of the results are to be Jde with conwJenL 


^•^•3 Speed of Computation 

An estimation of the speed of computation for any modeling effort of a 
co^resbor is dependent upon several variables: the number oAoluroes the 

minimum size of the volume, the number of sectors specified for a pr^allel! 
compressor representation, and to a lesser extent upon the number of per^L 

^r^sriHoS and sector. In the case of 

trie JBb 13 model, the compressor and combustor were sneclfieH ..c^nc lo i 

and up to 12 rlrcunfarnntlal sactnrs. The nlnl^rvru^ la^gSl^arlf^d'”' 
the maximum allowable time step in accord with the well known rpr /r 
Frledtlths, and Lawy) crltarla.. In this case. mO'" sL 

straints, typical run times on the GE/Honeywell 6000 computer were 167 time 
steps per minute for an undistorted- inlet model and 108 Jime steps per minute 
for a two-sector parallel compressor analyses. Quasi-steady-state thro^^W 

SSI "" “ d^Pdhdlng on the speed line 


Cour.mt, R., Friedrichs, K.O., and Lewy, 
zenglelchiingen der Matlicmatischen Physik 
pp 12-74. 


H., "Uber die Partiellen Dlfferen- 
, Math. Ann., Vol. 100, 1928, 


In this classical work, it was observed that a necessary condition for the 
convergence of a difference scheme is that the rate of propagation of signals 

Signal s^eed""'" " " 
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5.0 CONCLUSIONS AND RECOMMENDATIONS 


Accurate simulation of the clean-inlet-flow speed lines requires proper 
stacking of the individual stages of a compressor. In some instances, this 
requires modifying to varying degrees the level and/or the slope of the charac- 
teristics obtained from interstage test data. This type of modification 
requires the experience and judgement of the compressor designer to establish 
credible characteristics. It became apparent in the course of stacking the 
stages of the J85-13 compressor to obtain the overall compressor characteris- 
tics, that the pressure ratio-corrected flow speed lines passed through a 
maximum at the experimentally determined clean inlet surge line when the model 
reproduced the test surge line. In any region where the speed line slope was 
negative, the model was inherently stable. When the speed line slope became 
positive , the model was inherently unstable . 

The stability criterion used throughout this effort was based upon the 
ratio of the time rate of change of flow within a blade row to the time rate 
of change of the flow at the exit imposed throttling boundary condition. Where 
this value exceeded two in all blade rows, the model could not recover to stable 
performance. In all cases of clean inlet flow for both the "Moss" and "Mehalic" 
engines, this criterion led to consistent results and the critical value 
occurred In a pressure ratio-corrected flow region where an experimental aero- 
dynamic instability was observed. Further, this criterion produced consistent 
results throughout the distorted flow studies. For these reasons, it is 
assumed that the stability criterion is potentially an accurate predictor of 
compressor aerodynamic instability in a one-dimensional model for slow tran- 
sients. Further studies using other compression components will be required 
to establish the generality of this stability criterion. 

The search for a stability criterion led to another interesting result. 
Examination of the ratio of the exit flow from a blade-row volume to the 
entrance flow to that blade-row volume for each volume gives an Indication of 
the location where the flow breakdown within the compressor is occurring. As 
3 function of time , this parameter would Indicate values less than one in the 
stages upstream of the stage where flow breakdown was occurring and a value 
greater than one downstream of the stage where flow breakdown was occurring. 

This character is associated with flow storage in a blade-row volume and flow 
emptying from a blade— row volume, respectively. 

A dynamic blade row parallel compressor model with arbitrary extent 
sectors and with unspecified flow split was constructed. This model, using 
the clean inlet flow stability criterion applied to each sector, correctly 
predicted the trends of all the imposed 180", 1/rev distortions - total-pres- 
sure, total- temperature and combined total-pressure and total-temperatue. 

The accuracy of the distortion predictions in terms of pressure ratio/cor- 
rected fflow coordinate differences between predicted and measured values 
was quite good at high speeds where the speed lines were steep. At low cor- 
rected speeds where the speed lines had a much lower slope, the dynamic 
blade-row parallel-compressor model over-predicted the loss in surge pressure 





ratio. Classical steady-state parallel compressor theory is also known to 
overpredict the loss in surge pressure ratio when the speed line has a low 
value of slope. 

There are a number of factors to consider in discussing the differences 
between the distorted flow model predictions and the test data. First, the 
gap-to-radius ratios for the J85-13 compressor are quite small and thus, it is 
hard to imagine that significant inter-blade row flow redistribution is taking 
place. Secondly, since significant differences in the amount of rotation a 
streamtube associated with each sector experiences do not occur, it is assumed 
that the thermodynamics and velocity triangles are not major contributors to 
flow redistribution within the blade rows. 

As previously discussed, probably the most important source contributing 
to this difference lies in the unrealistic boundary conditions which are 
Imposed upon the flow at the compressor IGy. The analytical results of other 
Investigations clearly indicate that the pumping characteristics of the com- 
pressor can substantially modify the flow field between the measurement plane 
and the engine face by reducing the flow differences in the sectors, by estab- 
lishing a static-pressure distortion, and by establishing a swirl velocity 
component. Although static pressure differences between sectors (indicative of 
static pressure distortion) will be developed in this model upstream of the 
IGV s due to the flow split, none of the other effects are currently accounted 
for in the volumes between the measurement plane and the IGV's. Another possi- 
ble source contributing to the difference between predicted and measured loss 
in surge pressure ratio may lie in the fact that the third, fourth, and fifth 
stage bleeds are manifolded in a manner that would allow bleed recirculation 
from the high-static-pressure distorted sector to the low-static-pressure 
distorted sector. This possibility has not been taken into account during 
this program, although a very crude and preliminary approximation indicates 
that it is probably a secondary effect. Further analysis is required to 
fully determine the magnitude of the effect of bleed flow manifold recircula- 
tion on the loss in surge pressure ratio. 

The Dynamic Digital Blade Row Compression Component Stability Model is 
proving Itself as a viable tool for use in stability studies because of its 
to provide insight into the dynamic events at the surge line, the 
effects of circumferential distortion, and the stage where an aerodynamic 
instability originates. It also provides the detailed row-by-row variations 
in the state properties and vector diagram conditions throughout the com- 
pressor as they vary with engine operating conditions and inlet distortion. 
Furthermore, this approach to distorted inlet modeling has the potential to 
demonstrate whether the circumferential distortion sensitivity of a compression 
system can be determined analytically from clean inlet stage data plus a mini- 
mal amount of distorted inlet testing to verify the validity of the model. 

This procedure would be in contrast to the current practice of running exten- 
sive distorted inlet tests, followed by extensive development and manipulation 
of distortion parameters to correlate the resulting loss in surge pressure 
ratios. However, its main strength lies in the general form and flexibility 
of the equations to handle more complicated flow situations without having to 
resort to a new model or unduly alter the present model. 
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For these reasons, it is recommended that the present Dynamic Digital 
Blade Row Compression Component Stability Model be modified to include tan- 
gential flow redistribution in a gross sense, while maintaining the essence 
of the parallel compressor concept. In this manner, the compressor will 
establish more realistic inlet flow conditions at the IGV as opposed to the 
artificially established boundary conditions imposed at the measurement 
plane and assumed to be viilid back to the IGV. In addition, the effects of 
bleed flow redistribution within the bleed malnfolds should be investigated 
to determine the magnitude of this effect. 

The present model lacks the ability to predict stator induced instabil- 
ities since stator loss coefficients are assumed to be zero and the deviation 
angles are constant. Important insight could be obtained if realistic stator 
loss coefficients were determined and the stages restacked to give new rotor 
loss coefficients and deviation angles. Clean- inlet -speed-line and distorted- 
speed-line throttling would provide results which could be compared with the 
results of this study, specifically as to the location of the stage where the 
instability is Initiated and to determine if there are any differences in the 
manner in which flow breakdown occurs. 

The ability of the model to handle distortions with extents different 
than 180 “ is partially handled by allowing for flow redistribution upstream 
of the IGV’s but would also require that a method be developed for Including 
unsteady flow effects in the loss coefficients and deviation angles, rather 
than assuming that the loss coefficients and deviation angles determined 
from steady-state data are always applicable. 
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APPENDIX A 


STAGE CHARACTERISTICS 


This appendix documents the compressor stage characteristics 

utilized for this program. Characteristics for both the ’’Moss” and ’'Mehalic" 
versions of the engine are provided. Figures 45 through 58 illustrate the 
nondimens ional characteristics in terms of work and pressure coefficients as 
a function of flow coefficients. It should be noted that the nondlmens ional 
characteristics presented herein are normalized by pitch line wheel speed. 

The corresponding loss coefficient and deviation angle characteristics are 
presented in Figures 59 through 72, Tables 9 through 15 provide a tabulation 
of the coefficients used in the fourth-order polynomial representation of the 
loss coefficient and deviation angle characteristics. In addition, the bound- 
ing values of tangent of incidence angle over which the polynomial represen- 
tation is valid are indicated. 
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Figure 46, Preliminary and Final Pressure Coefficients, "Moss" Engine 
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Figure 50. Preliminary and Final Pressure Coefficients, "Moss" Engine 
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Figure 51. Preli-ii^ry and Final Work Coefficients, "Moss" Engine 
lOO^t N/V0. 
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Figure 52 . 
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Figure 58. Preliminary and Final Pressure Coefficients, »Mehalic" 
Engine 10096 N/v®. 
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Figure 65. Rotor Total -Pressure Loss Coefficients, "Moss" Engine 

1009 i N/v^. 
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Table 9. Polynomial Representation of Characteristics, "Moss" Engine, 
80% N/\/ 0 . 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, = Pn(TAN(i)) 

m' = Cj^TAN(i)^ + C2TAN(i)^ C3TAN(i)^ C4TANU) -f Cg 



X — 

TAN(i) DOMAIN 

POLYNOMIAL COEFFICIENTS 


ROTOR 

MAJC 

MIN 

Cl 

Ca 

C3 

n 


1 

.125 

.05 

10296.47 

-3725.04 


-32.3017 


2 


.08 


-1777.74 

375.2066 

- 34.4046 

1.242685 

3 


.02 

2474.917 

-649.511 

68.53722 

-3.14094 


4 

-095 

.025 


-345.946 


- 4.23584 


5 

-.005 

-.085 


-143.067 

-4.92062 

.6828110 


6 

.04 

-.10 

0.0 

-1.40873 

4.800277 

.5172106 


7 

.07 

-.10 

132.8065 

13.14993 

2.210511 

.2179109 

. 0664 l 4 l 6 

8 

.15 

0.0 

186.6458 

-20.3605 

-4.93^6 

.9791339 

.0078742 


B) TANGENT OF DEVIATION ANGLE, TAN( 6 ) = P^(XAN(i)) 


TAN(6) = D^TAN(i)^ + DalANU)^ D3TAN(i)^ D4TAN(i) ^ Dg 



TAN(i) DOMAIN 

POLYNOMIAL COEFFICIENTS 


ROTOR 

M/OC 

MIN 


1^2 



D3 

^4 

°5 

1 

. 125 

.05 

illRISlII 





2 

. 15 S 

.085 






3 

.09 

0.0 

2710.453 

- 446.232 

20.97019 

1.336071 

.0732273 



4 

.10 

0.0 

1101.343 

-163.322 

12 . 12245 

.5626970 

.1312668 

5 

-.005 

-.085 

8019.84 

1577.664 

110.2076 

2.964105 

. 1708325 

6 

.04 

-.09 

" 

0.0 

11 . 78 o 4 i 

4.899814 

WSSBSM 

.0690769 

7 

.07 

-.09 

555-3677 

15.65694 

.3223334 

.3663978 

.0978986 

8 

.14 

0.0 

-519.723 

163.2451 

-7.17728 


. 1256988 
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Table 10. Polynomial Representation of Characteristics, 
"Moss" Engine, N/v^ . 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, «»' = Fh(TAN(i)) 


U»' = C^TAN(i)^ + C2XAN(i)^ + C3TAN(i)^ + C4XANU) + C5 


ROTOR 

TAN(i) DOMAIN 

POLYNOMIAL COEFFICIENTS 

yiAX 

MIN 

^1 

^2 

C3 

^4 

^5 

1 

. 12 

.055 

9915.719 


569.3342 

-36.7901 

.9491094 

2 

.1325 

.06 

0.0 


64 . 16281 

-6.78050 

.3214183 

3 

.10 

• 02 

0.0 


23.51807 

-1.77353 

.1344891 

4 

.105 

.03 

0.0 

- 64.463 

27.43987 

-2.62657 

. 1695358 

5 

.005 

-.075 

0.0 

38. 15494 

16.47705 

1.266251 

. 1232226 

6 

• 04 

-.09 

0.0 

0.0 

7.049729 

.5873602 

.0974238 

7 

.07 

-.08 

0.0 

0.0 

3.781624 

.3201148 

1 ^ 

.0629112 

8 

.16 

-.02 

481.2131 

- 141.206 

12.46208 

.0541908 

.0224596 


B) TANGENT OF DEVIATION ANGLE, TAN( 6 ) = P^(XAN(i)) 


TAN( 6 ) = DjTAN(i)^ + DgTANU)^ + D^TAN(i )2 + D^TANU) + 


Dr\TV\D 

TAN(i) 1 

DOMAIN 

POLYNOMIAL COEFFICIENTS 

KUlv/K 

— 

MAX 

MIN 

•’l 

1 

‘>2 

“3 

*>4 

^5 

1 


.04 

14032.75 

-4379.77 

499.8880 


.5066225 

2 


.04 


-2314.34 

296.4176 


.4250809 

3 

. 11 

0.0 

6220.583 

-1438.49 

120.0122 


. 1375175 

4 

.12 

0.0 

2035.029 

-544. 69 

53.15900 


. 1441326 

5 

.02 

-.09 

2446.596 

CM 

00 

. 

00 

12.01995 

. 1958777 

.I 5040 l 8 

6 

.04 

-.075 

2617.987 

198.7423 

-.868978 


.0918752 

7 

.09 

-.10 

0.0 

0.0 

.7926447 

.134846 

.1157841 

8 

. 16 

0.0 

0.0 

-7.57951 

9.093984 


. 139405 


119 























































































































Table 11 


. Polynomial Representation of Characteristics, 
"Moss" Engine, N/\^. 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, I*) - Pn(lAN(i)) 


■JJj' _ c^TAN(i)^ + C2'lAN(i)^ C3lAN(i)^ C4lAN(i) C^ 


R0TX» 

1 ' ^ 

TAN(i) DOJttlN 

POLYNOMIAL COEFFK 

:iENTS 


M/OC 

MIN 

Cl 

C2 

^91 

mm 

^5 

1 

.13 

.04 

0.0 

0.0 

HIBSI 


. 1081784 

2 

.17 

.02 

0.0 





3 

.115 

.02 

0.0 





4 

.10 


0.0 





e 

-.0125 


0.0 

BBMM 

-11.9024 


. 1739506 

J 

6 

0.0 


10888.61 

-113.910 

-36.4737 


. 1360617 

7 

•o4 

-.07 

0.0 

-39.0041 

3.952665 


.0912396 

—L 

8 

.13 

.04 

0.0 

-6.74316 

2.590046 

.6832056 

-.005907 


B) TANGENT OF DEVIATION ANGLE, XAN(5) - Pj^(TAN(i)) 


TAN(6) = D^TAN(i)^ + D^TANU)^ D3TAN(i)^ D^TANCi) ^ D3 



^ ' X 

TAN(i) DOMAIN 

POLYNOMIAL COEFFK 

:iENTS 


ROTCn 

M/OC 

MIN 


^^2 

D3 

^4 


i 

.l45 

.05 


IBSSSm 

159.2616 



0 

.16 

.03 




WSS^M 

.3751164 

lb 

3 

.12 

.03 

0.0 

0.0 


. 1606487 

. 1324595 

4 





1630.419 


1.981923 



-.04 


284.6682 

42.44993 



A 

— 

.005 

-.03 

-24804.7 

655.0748 

124. 1600 

. 8608806 

.0921745 

7 

.05 

-.05 

0.0 

0.0 


.845641 


^ - 

8 

.13 

.06 

-2383.20 


-137.480 

7.178056 

1 .0309516 1 
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Table 12. Polynomial Representation of Characteristics, 
"Moss" Engine, 100% N/i/§. 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, «)' _ p„(TAN(i)) 

U)’ = Cj,TAN(i)^ C 2 lAN(i )3 CjTANU)^ C4lAN(i) -t- C5 



TAN(i) 

DOMAIN 

POLYNOMIAL COEFFICIENTS 

ROTOR 

M^X 

MIN 

— - 

Cl 

C2 

C3 

0 


— 1 

• 06 

.04 

0.0 

0.0 

215.8172 



2 

.082 

.056 

0.0 

0.0 




3 


.025 

0.0 

0.0 




4 

•11 


0.0 

0.0 

iBBHi 

mm 


5 

.02 



396.0822 

imm 



6 

.025 

• 

0 

32449.12 

313.5944 

-15.4741 


. 1240465 

7 

.09 

-.04 



2.912836 


. 1027257 

8 

.16 

.07 

5484.248 

-2612.16 1 

455.7857 


1.015377 


B) TANGENT OF DEVIATION ANGLE, TAN(6) = P^(TAN(i)) 


TAN(6) = D^TAN(i)^ + DglANCi)^ , D^IANU)^ + D^lAN(i) + 


DmV\D 

TAN(i) 

DOMAIN 

- ^ 

POLYNOMIAL COEFFICIENTS 

riAjxvut 


MIN 

Dl 


D3 

°4 

°5 

1 1 

.06 


Ol <2 

0.0 




2 








3 






bbb^b 


4 

WtSBM 






2.0016S1 

3 

0.0 

IbBSl 




1.850652 

■lim 

6 

.02 

-.01 

0.0 

144.509 


.0791817 



. 1128'527 

7 

.07 

-.06 

0.0 

0.0 




8 

.16 

.05 

0.0 

0.0 

13.39909 


.2463691 
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T^ble 13. Polynomial Representatiori of Characteristics, 
"Mehalic" Engine, 87 % N/v®. 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, U> = Pn(TAN(i)) 


U5’ = C^IAN(i)^ ^ CalANCi)^ ^ CjTANCi)^ C4TANti) C3 



X ^ 

TAN(i) DOMAIN 

POLYNOMIAL COEFFICIENTS 


ROTOR 

MWC 

MIN 


^2 

C3 



1 

.16 

.07 




BHH9 


2 

.22 

.05 

0-0 





3 

. 11 

.015 


warn 




4 

. 10 

.01 






5 

0.0 

-.10 





WS^^Si 

6 

.02 

-.12 

0.0 



Bi^SS 


7 

.05 

msm 

0.0 

RIRiSi 

1 6.686652 


.0850481 

8 


mm 

2.082061 

.2172444 



.0486026 


B) TANGENT OF DEVIATKMJ ANGLE, TAN(6) = Pj^(TAN(i)) 


TAN(6) = D^TAN(i)^ + DaTANCi)^ ^ D^TANCi)^ ^ D4TAN(i) ^ Dg 



TAN(i) DOMAIN 

POLYNOMIAL COEFFICIENTS 


ROTOR 

M4X 

MIN 

Dl 

»2 

D3 

°4 

°5 

1 

.16 

.09 

EBBS 

133-031^ 

-49.3825 

6.054828 

-. 172610 

0 

.22 

.065 

864. 1679 

BjBBi 


-2.59289 

.115214 

3 

.13 

-.04 

910.8485 

-190.641 

9.851025 


. 1504644 

4 

.09 

-.01 

0.0 

73.4386 

-12.8846 

.9956342 

.0577892 


.01 

-.075 

0.0 

94.01799 

14.72459 

mCSS^ 


^ 

6 

.02 

-.08 

3918.553 

545.4624 

21.88857 

lEiBSS 

.0340064 

7 

.045 

-.075 

2056.583 

103.7007 

-2.00420 

.1907098 

.0981537 

8 

.13 

-.02 

0.0 

2.627075 

-3.21001 


.0720412 
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Table l 4 . Polynomial Representation of Characteristics 
"Mehalic" Engine, 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, U)' = Pn(TAN(i)) 


m’ = C^TAN(i)^ -f C2TAN(i)3 CjTANd)^ C4mN(i) C.^ 



TAN(i) 

DOMAIN 

POLYNOMIAL COEFFICIENTS 

ROTOR 

M/UC 

MIN 

^1 

^2 

C3 

^4 

^5 

1 

.12 

.035 

0.0 

-50.1735 

21.89682 

-1.49564 

.1020144 

2 

.125 

.03 

664 . 1875 

- 353.044 

64.25218 

- 4.38497 

.1993741 

3 

.09 

-.02 

2740 . l 46 

-407.436 

27.85772 

-.634222 

.1037031 

4 

.07 

-.02 

4304.548 

-542.209 

33.16342 

-.242246 

.1166197 

5 

.02 

-.12 

1825-955 

374.4205 

34.43193 

2.633864 

.2201589 

6 

0.0 

-.12 

0.0 

0.0 

7.748663 

1.459335 

. l 64 o 66 o 

7 

.03 

-.08 


125.5993 

6.367800 

.7229095 

.099^66 

8 

.15 

-.02 

795.0528 

-241.583 

27.60512 

-.675066 

.0675536 


B) TANGENT OF DEVIATION ANGLE, TAN( 6 ) =P^(XAN(i)) 


TAN(6) s D^TAN(i)^ + DgTANCi)^ + DjTAN(i)^ + D^TANd) + 



TAN(i) 

DOM/VIN 

POLYNOMIAL COEFFICIENTS 


MAX 

MIN 

Dl 

»2 

D3 

“4 

^5 

1 

.12 

.05 

lISIEBfl 



-.945612 

. 1321859 

2 

.125 

.05 

0.0 

115.2502 


2.235098 

.0375340 

3 

.09 

-.04 

1288. 390 

-158.666 

6.088931 

1.075653 

. 1479905 

4 

.065 

0.0 

3717.266 

-478.975 

20.85731 

-.305552 

.0946924 

g. 

0.0 

-.09 

0.0 

90.67292 

23.07326 

1.358319 

.0974262 

6 

.02 

-.12 

33.11539 

5.302677 

-.411678 

-.193495 

.0429261 

z 

.07 

-.10 

720.2954 

31.11858 

-1.48927 

.0869565 

.0774906 

8 

.15 

-.03 

358.2427 

.♦70.2615 

4.949831 

-.146842 

.0644381 
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Table 15. Polynomial Representation of Characteristics, 
"Mehalic" Engine, 1CX)?6 N/v 9 . 


A) RELATIVE TOTAL-PRESSURE LOSS COEFFICIENT, «>' = Pn(TAN(i)) 


IS’ = C.,TAN(i)^ ^ CalAN(i)^ t C3TAN(i)^ ^ C4lAN(i) ^ Cg 


ROTOR 

1 T 

TAN(i) DOMAIN 

POLYNOMIAL COEFFK 

:iENTS 


M(OC 

MIN 

Cl 

C2 

C3 

mm 


1 

.06 

.02 

0.0 





2 

-o 4 

-.01 

0.0 



HRSIB 


3 

.03 

-.05 

0.0 





4 

-OS 

-.03 

0.0 






-.06 

-.l 4 

0.0 





6 

-.01 

-.13 

0.0 

0.0 

2«Q2 u 1 l 1 sIsasI 

1.514769 

. 2026557 

7 

•04 

-.08 

0.0 




. 1210497 

8 

.13 

0.0 

0.0 

1 0.0 

16 . 105443 


. 1204295 


B) TANGENT OF DEVIATION ANGLE, TAN(6) = Pj,(TAN(i)) 


TAN(6) = DiTAN(i)^ ^ D2lAN(i)3 > D3TAN(i)^ ^ D^TANU) > D3 


ROTOR 

— ^ - j, 

TAN(i) DOMAIN 

POLYNOMIAL COEFFK 

:IENTS 


ig 

MIN 

Dl 

»2 

D3 

^4 

•^5 

1 

URS^H 

.025 


-l46o6.7 



.4784627 

p 

•05 


66410.15 

-3616.78 




1 

.09 

-.08 


-4.16092 

-6.00985 

-.019592 

.1019519 

4 

.075 

-.06 

2586.269 

-134. 166 

-7.2l4l8 

.9417444 

.0524781 


0.0 

-.1 

0.0 

134.3838 

28.72648 

1.115973 

.0368576 

^ 

1 

.04 

-.09 

0.0 

19.85554 

2.869345 

-.115185 

.0048145 

u 

7 

.05 

-.08 

0.0 

-1.69849 

. 8652017 

.3387322 

.0648323 

^ 

8 

.13 

-.01 

0.0 

16.42562 

-2. 87861 

.3745065 

.0179125 
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APPENDIX B 


LIFT DIRECTION CORRECTIONS ANGLES 


The figures and tables contained in this appendix document the lift- 
direction correction angles (gj.) used in this program. The distribution of 
the tangents of the correction angles for the rotors and stators along a 
speed line has been represented as a linear function of the tangent of the 
incidence angle. Figures 73 through 79 Illustrate the final correction angle 
distribution as used in the dynamic model. Tables 16 throu^ 22 provide a 
tabulation of the coefficients used in the straight-line representation of 
the correction angles. The inlet guide vane and outlet guide vanes were 
specified in the modeling effort as having constant incidence angle over a 
speed line and as such prevent presentation of the lift direction correction 
angle as a function of incidence angle. Therefore, the lift direction correc- 
tion angles for the inlet and outlet guide vanes were chosen to be constants 
equal to the average of the calculated distribution of values. Table 23 
documents the values used for both engines at each corrected speed. 
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Figure 73 - Li 










KEPKOT)UC!Bn.!TY OP THE 
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Table l 6 . Lift Direction Correction Angle Coefficients, 
"Moss" Engine, 8 o 56 N/v 0 . 


TAN(Pc) = (M)(TAM(i)) + B 


BLADE 

ROW 

B 

(y-intercept) 

M 

(SLOPE) 

R1 

.016271 

.0047059 

SI 

-.019292 

-.200995 

R 2 

.022483 


S 2 

->037873 

-.19494 

R 3 

.024924 

-.066187 



-.047869 

-.1375 

R 4 

.018279 

-.0085714 

S 4 

-.041828 

-.14857 

R 5 

.019365 

• 1554 

S 3 

-.03591 

-.16696 

R 6 

i 

.010067 

.10489 

S 6 

-.03382 

-.2061 

B 7 

-.0008977 

.026047 

S 7 

-.029025 

-.14517 

R 8 

-.0019556 

-.0044444 

S8 

-.0427 

-. 15269 
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Xable 17. 













BLADE 

ROW 


B 

(Y-INTERCEPT) 


N 

(SLOPE) 





-.0005394 

-.04536 


.078571 

-.25607 

-.15484 

-.41935 

-.20333 

-.104 

-.19813 

-.12364 

-.074074 

-.21778 

.013913 

-.19434 

-.030252 

-.1587 

-.064789 

-.1657 
























Table 19- 














Table 20. Lift Direction Correction Angle Coefficients, 
"Mehalic" Engine, 8 ?% N/v^. 


lAN(Pc) = (M)(lAN(i)) + B 


BLADE 

ROW 

B 

( Y-INTERCEPT) 

M 

(SLOPE) 

R1 

.030918 

-.10256 

SI 

—Oiq 775 

-.21149 

R2 

.033216 

-.12327 

S 2 

-.037727 

-.16881 

R 3 

.020568 

—095385 


-.051514 

-.2015 

R 4 

.025781 

-.07205 

S 4 

-.047627 

-.24545 

R 5 

.0081652 

-.093939 



-.036741 

-.21875 

R 6 

.0024279 

-.018605 

S6 

-.034042 

-•21953 

R 7 

-.0055562 

-.071006 

57 

-.028844 

-.14444 

R 8 

-.0029805 

-.04878 

S 8 

-.04808 

-.22027 























































T^ble 22. Lift Direction Correction Angle Coefficients, 
"Mehalic" Engine, 1CX)96 N/v®. 


TAN(Pc) = (M)(TAN{i)) + B 


BLADE 

ROW 

B 

(Y -INTERCEPT) 

M 

(SLOPE) 

R1 

.01337 

-.074074 

SI 

-.014335 

-.084472 

R2 

.0152 

-.164 

52 

-.04528 

-.21 

R3 

.01123 

-.036667 

S3 

-.057699 

-.21078 

R4 

.015044 

-.055556 

54 

-.051818 

-.19273 

R5 

.012426 

.067532 

S5 

-.040244 

-22095 

r6 

-.0036409 

-.030108 

56 

- .036708 

-.21695 

R7 


-.034783 

57 

-.030629 

-.18295 

r8 

_.00712Q5 

-.039306 

S8 

-.049386 

-.2656 



















Table 23. Average tan (P^.) for IGV and OGV for "Moss" and 
"Mehalic" Engines . 


TAN(Pq) = (M)(XAN(i)) + B 


SPEED 

BLADE 

B 

M 

ROW 

(Y-INTEHCEPT) 

(SLOPE) 

MOSS 80% 

IGV 

.000190 

0.0 

MOSS 8096 N/^/§" 

OGV 

-.000645 

0.0 

MOSS 87 % nZ-T^" 

IGV 

.014013 

0.0 

MOSS 87 % N/^/§’ 

OGV 

-.000648 

0.0 

MOSS 9^ N/v/Sr 

IGV 

.283881 

0.0 

MOSS 9^!^ N/7^ 

•OGV 

-.000657 

0.0 

MOSS 100% N/^/§' 

IGV 

0.0 

0.0 

MOSS 100% N/^/e 

OGV 

-.000646 

0.0 





MEHALIC 87 % N/v^ 

IGV 

.040916 

0.0 

MEHALIC 87 % NZ-Z^ 

OGV 

-.060647 

0.0 

MEHALIC 94 % NZ>/0 

IGV 

0.0 

0.0 

MEHALIC 94 % NZv^ 

OGV 

-.000652 

0.0 

MEHALIC 100% NZZ^ 

IGV 

0.0 

0.0 

MEHALIC 100 % NZV^ 

OGV 

-.000647 

0.0 
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APPENDIX C 


CLEAN INLET DOCUMENTATION 


The computer output listing including herein serves to illustrate the 
type of documentation available for the clean inlet modeling of both the 
"Moss** and **Mehalic** J85-13 engines. In order to aid in the interpretation 
of the output, an explanation of the parameter titles is presented in Table 
24. It should be noted that output is provided only for volumes occupied 
by either stationary or rotating blade rows. 

Table 25 provides a tabulation of the compressor performance for the 
high flow condition on the "Moss" engine 100 percent speed line. 


141 



Table 24. Computer Listing Output Parameters. 


PCTNC 

WCORR 

P/P-OA 

EXIT FF 

CZl 

CUl 

WUl 

CT 

vrr 

U1 

M-ABS 

M-SEL 

ALPHA 

BETA 

PSI 

PTl 

TSl 

TTl 

TNI 

INC 

LOSS 

TND 

DEV 

DFACT 

PHI 

PSI 

PSI-P 

PRl 

PR2 

TRl 

TR2 


- % Corrected Speed 

- Corrected Inlet Flow (kg/sec) 

Overall Pressure Ratio 

- Exit Flow Function 

- Inlet Axial Velcoity (m/sec) 

- Inlet Absolute Swirl (m/sec) 

- Inlet Relative Tangential Velocity (m/sec) 

- Inlet Absolute Velocity (m/sec) 

- Inlet Relative Velocity (m/sec) 

- Inlet Pitchline Wheel Speed (m/sec) 

- Inlet Absolute Mach Number 

- Inlet Relative Mach Number 

- Inlet Absolute Air Angle (deg.) 

- Inlet Relative Air Angle (deg.) 

- Inlet Static Pressure (N/cm^) 

- Inlet Total Pressure (N/cm^) 

- Inlet Static Temperature (“ K) 

- Inlet Total Temperature (® K) 

- Tangent of Incidence Angle 

- Incidence Angle (deg.) 

- Total-Pressure Loss Coefficient 

- Tan g ent of Deviation Angle 

- Deviation Angle (deg.) 

- Diffusion Factor 

- Flow Coefficient 

- Work Coefficient 

- Pressure Coefficient 

- Cumulative Pressure Ratio 

- Blade-Row Pressure Ratio 

- Cumulative Temperature Ratio 

- Blade-Row Temperature Ratio 
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Table 24* Computer Listing Output Parameters (Concluded) 


AD-EF 

W2/W1 

DWX 

DW/DWEX 


Stage Adiabatic Efficiency 

Ratio of Exit Flow to Inlet Flow 

Volume Averaged Flow Time Derivative (kg/sec^) 

Ratio of Volume Averaged Flow Time Derivative to Exit 
Volume Averaged Flow Time Derivative 



Table 25* Calculated Performance, Clean Inlet Flow, "Moss'» Engine 100% N/n/6 
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APPENDIX D 


TOTAL-PRESSURE DISTORTION DOCUMENTATION 


This appendix contains information which illustrates the type of docu- 
mentation available for the 180°, 1/rev distortions modeled for the **Moss” 
J85-13 engine. Figues 80 through 82 present the distortion profiles as 
deduced from the distortion instrumentation measurements for three levels 
of distortion at 100 percent corrected speed. Both radial and circumferential 
profiles are supplied. The normalizing average pressures are indicated as 
well as the NASA reading number from which the data was obtained. Reference 
4 contains a complete description of the number, type, and location of the 
instrumentation probes. 

Detailed compressor maps of the total-pressure distortion, dynamic 
parallel-compressor analyses are presented in Figures 83 through 85. Included 
on tiu* figures is the performance of the individual parallel-compressor sectors 
and the resultant overall performance. 

Documentation of the distortion cases analyzed is presented in the form 
of computer output listings for which the reader is referred to Table 24 
and supplemental Table 26 for an explanation of the parameter titles. A tab- 
ulation of compressor performance for each of the sectors for the high level 
total-pressure distortion at 100 percent speed is supplied in Table 27. 
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IMMERSION 5 




1.0 
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ANGULAR LOCATION, AFT LOOKING FORWARD, DEGREES 
(0® = TOP DEAD CENTER) 



(a) Circumferential Profiles. 

Figure 80. Circumferential Total-Pressure Distortion Profiles 
(RDG 478 ), ’*Moss" Engine 100% N/v^. 

















IMMERSION 5 


O o O O 


IMMERSION 4 


o O O o o o 


IMMERSION 3 


_[o 0^000 


IMMERSION 2 




0 

0 

0 

0 

0 

0 0 0 0 0 0 



IMMERSION 1 

O O 


O O 


Jo o o.Q o O. ^ . 
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ANGULAR LOCATION, AFT LOOKING FORWARD, DEGREES 
(0° = TC»> DEAD CENTER) 


(a) Circumferential Profiles. 

Figure 82 . Circumferential Total -Pressure Distortion Profiles 
(RDG 99), "Moss” Engine 100% 
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Figure 82 . 
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CLEAN INLET SURGE LINE 
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CORRECTED FLOW (kg /sec) 

Figure 85 . Effect of l80°, l/Rev Total -Pressure Distortion (9M Screen) on Surge Line of 
”Moss” Engine - Sector Performance, 


Table 26. Supplemental List of Con^uter Output Parameters. 


PCTNCL 

ANGLE 

TALCW 

P/P-OA(L) 

T/T-OA(L) 

DPTR-IN 

DTTR-IN 

DPSR-IN 

ROT 

DPTR2 

DTTR2 

DPSR2 


Sector Local Corrected Speed 
Angular Extent of Sector 

Equivalent Total Area Sector Corrected Flow (kg/sec) 
Sector Overall Pressure Ratio 
Sector Overall Temperature Ratio 

Normalized Total-Pressure Anqilitude at Inlet, iSPj/Pj 
Normalized Total-Temperature Amplitude at Inlet, AT^p/Tj 
Normalized Static-Pressure Amplitude at Inlet, AP 5 /FS 
Cumulative Sector Rotation (deg.) 

Blade-Row Exit Normalized Total-Pressure Amplitude 
Blade-Row Exit Normalized Total-Temperature Amplitude 
Blade-Row Exit Normalized Static -Pressure Amplitude 
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APPENDIX E 


total-temperature distortion documentation 


In the same manner as Appendix D, this appendix illustrates the detailed 
supplemental information that is available as model output. Figure 86 pre- 
sents the circumferential and radial profiles of the temperature distortion 
as indicated by the NASA test data. Included on the figure are the NASA 
reading number and the normalizing parameters. A detailed compressor map of 
the temperature-distortion throttling simulations, illustrating the oP^atlng 
points of the parallel-compressor sectors is presented in Figure 87 Table 
28 presents documentation of the compressor performance for the -Mehalic 
engine 100 percent corrected flow high-flow condition. 
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IMMERSION 5 



IMMERSION 3 




IMMERSION 1 



ANGULAR LOCATION, AFT LOOKING FORWARD, DEGREES 
(0° = TOP DEAD CENTER) 


(a) Circumferential Profiles. 

Figure 86. Circumferential Total-Temperature Distortion Profiles 
(RDG 568), "Mehalic” Engine 86.8% N/^/w, 
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CORRECTED FLOW (kg /sec) 

Figure 8?. Effect of l80°, l/Rev Total- Temperature Distortion on Surge Line of "Mehalic" 
Engine - Sector Performance. 
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COMBINED TOTAL-PRESSURE AND TOTAL-TEMPERATURE DISTORTION DOCUMENTATION 


Sample radial and circumferential profiles of the 180°, 1/rev total- 
pressure and total-temperature distortions as taken from the NASA data for the 
opposed, coincident, and 90° overlapped orientations are presented in Figures 
88 through 90. The normalizing parameters and NASA test reading numbers are 
indicated on the plots. 

Figures 91 and 92 represent the opposed orientation. Figure 93 the coin- 
cident orientation, and Figures 94 and 95 represent the 90° overlapped orien- 
tation with the operating points of each of the parallel compressor sectors 
shown on each figure. The local corrected speeds of each sector are also 
indicated on the maps. 

Table 29 illustrates the type of documentation that is available as model 
output, in this case, for the "Mehalic" engine 100 percent corrected flow 
condition. 
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IMMERSION 5 



IMMERSION 3 



IMMERSION 2 



IfWERSION 1 



ANGULAR LOCATION, AFT LOOKING FORWARD, DEGREES 
(O® = TOP DEAD CENTER) 


Figure 88* 


(a) Circumferential Temperature Profiles. 

Combined Circumferential Total-Pressure and Total- 


Temperature Distortion Profiles, Opposed Orientation 
(RDG 479), "Mehalic" Engine 92 . 59^ N/v^. 
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(b) Radial Temperature Profiles. „ 

Figure 88. Combined Circumferential Total -Pressure and Total -Temperature Distortion Profiles 
Opposed Orientation (RDG 4?9), "Mehalic" Engine 92.5% N/V« (Continued). 
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(c) Pressure Profiles. 

Figure 88. Combined Circumferential Total-Pressure and Total- 

Temperature Distortion Profiles, Opposed Orientation 
(RDG 479), "Mehalic" Engine 92.5% N/v/S (Concluded). 
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Figure 89 • 


(a) Circumferential Temperature Profiles. 

Combined Circumferential Total-Pressure and Total- 
Temperature Distortion Profiles, Coincident Orientation 
(RDG 423), "Mehalic" Engine 99 * 15 ^ N/v . 
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(b) Radial Temperature Profiles. 

Figure 89 . Combined Circumferential Total-Pressure and Total-Temperature Distortion Profiles, 
Coincident Orientation (RDG 423), "Mehalic" Engine 99.1% N/\/9 (Continued). 
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(c) Pressure Profiles. 

Figure 89. Combined Circumferential Total-Pressure and Total- 

Temperature Distortion Profiles, Coir^ident Orientation 
(RDG 423), "Mehalic" Engine 99»1% N/v® (Concluded). 
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IMMERSION 5 



IMMERSION 2 

1.05p 

1.00 — Q — Q — Q. 

O O O o o o 

.95 - 


IMMERSION 1 



ANGULAR LOCATION, AFT LOOKING FORWARD, DEGREES 
(0° = TOP DEAD CENTER) 

(a) Circumferential Temperature Profiles. 

Figure 90 . Combined Circumferential Total -Pressure and Total- 

Temperature Distortion Profiles, 90^ Overlapped Orientation 
(RDG 420), "Mehalic” Engine 99.0% N/TQ. 


171 







ANGULAR LOCATI(»I, AFT LOOKING FORWARD, DEGREES 
(0° = TW DEAD CENTER) 


5 


t/3 

5 

M 

Q 




30® RAKE 
150° RAKE 

2100 rake 
3300 rake 


(c) Pressure Profiles. 

Figure 90. Combined Circumferential Total-Pressure and Ibtal- 

TMperature Distortion Profiles, 90® ^erlapped Orientation 
(RDG 420), "Mehalic" Engine 99.0% N/y/® (Concluded). 
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Fiflure 95 Effect of 18 qo, l/Rev Combined Total-Pressure and Total -Temperature Distortion 
90° Overlapped Orientation Moderate Total-Temperature Distortion on Surge Line 
of "Mehalic" Engine - Sector Performance. 
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(c) Sector 3« 

Table 29. Calculated Performance, Combined Total-Pressure and Total-Temperature 
Distortion (RDG 420), "Mehalic" Engine 99-0% N/v^ (Continued). 
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APPENDIX G 


DISTORTION TRANSMISSION DOCUMENTATION 


identifying regions of amplification and attenuation of 
distortion in the compressor, the distortion amplitudes at the exit of each 

AjSdL~s r 1„ the listlnge preset 

«=rK?7 appendix is a compilation of the amplification results 

a lished the distortion analysis performed at the intermediate flow 
Sen normalized distortion amplitudes are formulated as the dif- 

bv thr maximum and minimum values normalized 

^ota? nrlsr®^ parameter values. Plots are provided for normaUzed 

t tal pressure, total-temperature, and static-pressure amplitudes. 
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Figure 97. ture *«pll,ic.ti„„ For Inlet lot. 1-Preesure Dietortlon 
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Figure 99. Predicted Total-Pressure Amplification For Inlet Total-Temperature Distortion. 
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Figure 101. Predicted Static-Pressure Amplification For Inlet Total-Temperature Distortion. 
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Figure 104 . Predicted Static-Pressure Amplification for Combin 
Total -Temperature Distortion, Opposed Oreintation. 
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Figure 105. Predicted Total-Pressure Amplification for Combined Inlet Itotal-Pressure and 
Total-Temperature Distortion, Coincident Orientation. 
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Figure 105* Predicted Total— Pressure Amplification for Combined Inlet Total -Pres sure and 
Tota 1- Tempera tu re Distortion, Coincident Orientation. 
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Figure 106, Predicted Total -Temperature Amplification for Combined Inlet Total-Pressure 
and Total-Temperature Distortion, Coincident Orientation. 
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Figure 107. Predicted Static-Pressure Amplification for Combined Inlet Total-Pressure 
and Total-Temperature Distortion, Coincident Orientation. 
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Figure 109. Predicted Total-Temperature Amplification for Combined Inlet Total-Pressure 
and Total -Temperature Distortion, 90° Overlapped Orientation. 
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